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We systematically investigate the new experimental method of using atoms or molecules 
to measure important parameters of neutrinos still undetermined; the absolute mass 
scale, the mass hierarchy pattern (normal or inverted), the neutrino mass type (Majo- 
rana or Dirac), and the CP violating phases including Majorana phases. Most of these 
observables are difficult to measure in neutrino oscillation experiments. 

There are advantages of atomic targets, due to the closeness of available atomic ener- 
gies to anticipated neutrino masses, over nuclear target experiments such as the end 
point spectrum of (3 decay and two-electron line spectrum in the neutrinoless double /3 
decay, both of which address some of the overlapping objectives with atomic/molecular 
experiments. Disadvantage of using atomic targets, the smallness of rates, is overcome 
by the macro-coherent amplification mechanism. 

The atomic or molecular process we use is a cooperative deexcitation of a collective 
body of atoms in a metastable level |e) emitting a neutrino pair and a photon; |e) — >■ \g) + 
7 -|- UiVj where fj's are neutrino mass eigenstates. The macro-coherence is developed by 
trigger laser irradiation of two colors, which frequently causes two-photon process |e) -H- 
Iff) + 7 + 7 , |e) -|- 7 O |(?) -|- 7 inside the target. We discuss important aspects of the 
macro-coherence development in detail, by setting up the master equation for the target 
Bloch vector (whose components are population difference and medium polarization) and 
propagating electric field. Our master equation includes effects of phase decoherence of 
medium polarization and decay of population difi^erence. 

The spectral rate (the number of events per unit time) of macro-coherent radiative 
emission of neutrino pair has three parts, and is given by a factorized formula of the 
form, (overall oj independent rate denoted by Fq) x (spectral shape function denoted by 
/(cj)) X (time evolving dynamical factor), where lu is the photon energy. The constant 
factor Pq determines the overall rate in the unit of 1/time, and for Xe it is of order, 
1 Hz (n/10^^cm~^)^ (V/lO^cm^). The dynamical factor is time dependent and is given by 
the space integrated quantity over the entire target, of the product of magnitude squared 
of coherent polarization and field strength (in the unit of the maximally extractable 
energy density) stored inside the target. The asymptotic value of time evolving dynam- 
ical factor is given by contribution of field condensate accompanied by macroscopic 
coherence, which is calculated using the static limit of the master equation. With an 
appropriate choice of heavy target atoms or molecules such as Xe and I2 that has a 
large Mix El matrix element between |e) and \g), we show that one can determine 
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three neutrino masses along with distinction of the mass hierarcliy pattern (normal or 
inverted) by measuring the spectral shape I(uj). If one uses a target of available energy 
of a fraction of 1 eV, the most experimentally challenging observable, the Majorana 
CP phases, may be determined, comparing detected rate with differences of theoretical 
expectations which exist at the level of several %. The Majorana CP violating phase is 
expcctcid crucial to the understanding of the matter-antimatter imbalance of our uni- 
verse. Our master equation, when applied to ElxEl transition such as pH2 vibrational 
Xv = 1 — > 0, can describe explosive PSR events in which most of the energy stored in 
|e) is released in duration of order a few nano seconds. 

The present paper is intended to be self-contained explaining some details of related 
theoretical works in the past, and further reports new simulations and our ongoing 
experimental efforts of the project to realize the neutrino mass spectroscopy using 
atoms / molecules . 

Subject Index neutrino, Majorana particle, lepto-genesis, macro-coherence, paired super- 
radiance 
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1. Introduction and overview 

1.1. Remaining important problems in neutrino physics and our objective 

The present status of neutrino mass matrix is summarized by the following central values 

measured by oscillation experiments [1], [2]: 



sf2 = 0.31 , = 0.42 
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The usual notation of angle factors is used; Sij 
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for Majorana neutrinos , = 1 for Dirac neutrinos , (4) 



(5) 



(where Mo is the diagonalized mass matrix). Neutrino masses are ordered by > m2 > mi 
for the normal hierarchical mass pattern (NH) and m-2 > mi > m^ for the inverted hierarchy 
(IH). For convenience we define the smallest mass by mo, which is = mi for NH and = ms 
for IH. 

The ongoing and planned experiments to measure the neutrino masses using nuclei as 
targets are in two directions; (1) measurement of the beta spectrum near the end point 
sensitive to both Dirac and Majorana masses, (2) neutrinoless double beta decay near the 
end point of two electron energy sum, sensitive to Majorana masses alone. In the neutrinoless 
double beta decay one attempts to measure the following parameter combination called the 
effective neutrino mass [3]; 

2 
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^3*13 + ^2^12^13 + nT'iCu^^is + 2mim2Si2Ci2Ci3 cos(2a) 



+2mim3Si3Ci2Ci3 cos 2(^ — 5) + 2m2m3Si3Si2Ci3 cos 2(a — (3 + 5) 



(6) 



using our convention of Majorana phases. The best upper limit of neutrino mass scale is 
derived from cosmological arguments, and is 0.58cV (95% CL) [4]. 

Despite of this remarkable success in neutrino physics, there are still many important ques- 
tions to be answered. (1) Whether the nature favors either of the neutrino mass type, Dirac 
(described by 4 component spinor equation) or Majorana mass (described by 2 component 
spinor), is unknown despite of its vital importance to Icpto-gciicsis theory [5], [6]. (2) Two 
important parameters of the neutrino mass matrix Mu, the smallest mass and the CP 
violating phase 5, are inaccessible experimentally in the near future. (3) A definite principle 
of measuring the additional CP phases a, /3 [7], [3] intrinsic to the Majorana neutrino has 
not been proposed as yet. 

These are challenged by our method of neutrino mass spectroscopy [8], [9], [10] using atoms 
or molecules instead of nuclei as targets, as closely explained in the present article. 
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Fig. 1 A— type atomic level for RENP |e) — t- \g)+^ + i'iVj with vi a neutrino mass 
eigenstate. Dipole forbidden transition |e) — t- I^) + 7 + 7 may also occur via weak Mix El 
couplings to virtual intermediate state \p). 

In the rest of this Section we overview theoretical aspects by summarizing essence of 
theoretical sections, Sec. 2 and Sec. 3. Experimental status of our project is summarized in 



1.2. Radiative emission of neutrino pair (RENP) 

The atomic process we use for experiments is |e) — )• l^f) + ^ + uiUj, as illustrated in Fig.l. 
The process exists in the standard electroweak theory [11] without any doubt: it occurs 
as a combined effect of weak interaction and QED, a kind of second order perturbation 
process when one regards four Fermi weak process as the first order process. We denote 
this process by RENP (radiative emission of neutrino pair). The radiation- less pair emission 
|e) — 7- + ViVj is faster as an elementary process than RENP, but with the aid of macro- 
coherent amplification by trigger laser irradiation RENP dominates over the radiation-less 
process. Moreover, the single emitted photon is a key to obtain needed information on 
neutrinos, since emitted neutrinos are difficult to detect. The atomic state |e) is assumed 
metastable, which means to us that its lifetime roughly > 1msec (its optimal value to be 
determined by repetition cycle of excitation and trigger irradiation in actual experiments). 
Besides a single photon 7 the final state has two neutrino mass eigenstates of i^j , i = 1, 2, 3. 
The crucial key element of our experimental methods is the ability of resolving neutrino mass 
eigenstates rather than flavor eigenstates, as realized by the excellent frequency resolution of 
used trigger laser. This becomes possible by using the technique of trigger laser irradiation in 
atomic processes. There are step function like threshold rises in the photon energy spectrum 
and six threshold locations are at 



Here eah = £a — is the atomic energy difference between two states |a) , Threshold loca- 
tions sensitive to neutrino masses are separated by small photon energies. For example, 
{rrii + mj)^/(2eeg) ~ 5 meV for rrii + nij = O.leV and eeg = leV, taking a typical atomic 
energy difference. One can separate different mass eigenstates by fully exploiting the accu- 
racy of frequency in the range of oj < cou (the largest threshold), used as trigger laser in 
our proposed experiments. In our approach one does not need this order of precision of 
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detected photon energy. From the continuous spectral shape of the single photon rate one 
can determine both neutrino masses and mixing angles, along with Majorana phases a, (5. 

Atomic/molecular targets have an advantage over conventionally used nuclei in their close- 
ness of released energy to neutrino masses (expected much smaller than 1 eV) . A demerit of 
these targets is the weakness of RENP rate oc G^ae"^ with e^g the available energy of order 
eV and n ~ 5. The smallness of rate is due to the small Fermi constant Gp ^ 10~^^eV~^. We 
use macro-coherence (giving the rate v^V with n the target number density, V the target 
volume) to overcome this problem. When the number density n is close to the Avogadro 
number per cm^, RENP rate may become measurable. 

If the macro-coherent amplification works as expected, the neutrino pair emission accom- 
panied by the photon of energy oj occurs according to a time dependent rate formula of the 
factorized form, 

V^2u{uJ,t) = ToI{uj)r]^{t). (8) 

The constant RENP rate Fq may become of order 1 Hz at the target number density 
lO^^cm^^ (the Fq value scaling with the number density oc n^) and the volume lO^cm^ 
in the Xe example whose RENP spectrum I{oj) in the threshold region is shown in Fig. 2. 
^ Relevant levels of Xe are |e) = 5j3^(^P3/2)6s^[3/2]2, metastable with lifetime of O[40]sec, 
b) =^P^^Sq, and intermediate state \p) = 5p^(^P3/2)6s^[3/2]i. The initial and the final 
states, |e) , have different parities and the angular momentum difference AJ = 2. Excited 
states here may be described by a pair state of 6s electron and 5j9 hole: both of the initial 
metastable and intermediate states are given by the spin triplet pair with a large breaking 
of LS coupling scheme. The dynamical factor riui{t) shall be discussed later in Sec. 3 based 
on the solution of the master equation of Sec. 2. 

Prom the spectrum feature of Fig.2 it should not be difficult to measure the absolute 
neutrino mass scale and the distinction of normal (NH) and inverted (IH) mass hierarchy of 
neutrino mass pattern. The Majorana vs Dirac distinction is harder for Xe due to a large 
energy level difference e^g ~ 8.3 cV. It is found in [12] that the appropriate energy scale for 
measurements of Majorana CP phases is a fraction of eV. 

Attractive candidates of targets may be found in molecules. Molecules are interesting due 
to a rich vibrational and rotational band structure with much smaller level spacing than 
the electronic ones. We illustrate RENP spectrum calcTilation for I2 electronic transition in 
Fig.3. Difference of spectrum rates for different CP Majorana phases is a several to 10 % for 
IH. 

1.3. Paired super-radiance (PSR) to he controlled 

It is crucial for the success of our method to control a twin process, PSR (paired super- 
radiance), |e) — >■ Iff) + 7 + 7- PSR event is interesting in itself. With the macro-coherence 
(in which the coherent volume is not wavelength limited unlike the single photon super- 
radiance (SR) [13]), typically exhibiting the back to back two photon emission with equal 
energies at the half of atomic level difference eeg/2 (under the trigger laser of the same 
frequency). Emitted two photons are highly entangled especially for J = — >■ transition. 



^ Neutrino parameters taken for this calculation are somewhat different from the most recent values 
given above, Eq. 2 and are given in Sec. 3. 



5/85 



Xe, Dirac NH vs IH: mO=l,10,50meV 




Fig. 2 RENP dimensionless spectrum function /(w) near tlie neutrino pair emission 
thresholds from Xe level 5p^(^P3/2)6s^[3/2]2. Neutrinos of the smallest mass of 1, 10 and 
50 meV are taken for the normal (solid curve) and the inverted (dashed curve) hierarchical 
mass pattern. 



12 A'v=l->Xv=15: mO=20raeV 
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Fig. 3 I2 RENP spectrum between A' v=l and Xv=15. The Majorana (in colored solid 
for NH and in colored dashed for IH) vs Dirac (in black solid for NH and in black dashed 
for IH) cases are compared. Three Majorana CP phase combinations (a, f3 — 5) = (0, 0) (in 
red), (7r/2,0) (in green), and (0, 7r/2) (in blue) are taken, with the smallest neutrino mass 
20 meV. The vertical scale is in arbitrary units. 



For long targets of the number density of metastable state over ~ lO^^cm"'^, explosive PSR 
occurs if the initial coherence between |e) and \g) is present [14]. Explosive PSR event is 
characterized by instantaneous release of energy stored in the upper level |e) into short 
pulses of some time structure, taking place with a time delay after weak trigger irradiation, 
as illustrated in Fig. 4 for vibrational pH2 transition. The largest instantaneous output/input 
power ratio in this figure is O[10^^]. The macro-coherent PSR is fundamentally different from 
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Fig. 4 Time-evolving PSR output flux resulting from the symmetric trigger irradiation 
at two target ends, of the power range, 10^^^ ~ lWmm~^, under the conditions of the 
target number density n = 1 x lO^^cm"^, target length = 30cm, relaxation times T2 = T3 = 
10, Ti = 10^ ns's, and the initial polarization, rf^^ = 1 and all other Bloch vector components 
= 0. Depicted outputs from 1 Wmm~^ trigger power in black, from 10~^Wmm~^ in red, and 
from 10~^^Wmm~^ in blue are displaced almost equi-distantly in the first peak positions. 
Vibrational transition Xv = 1 — ?■ Xv = of pH2 is considered. 

SR of a single photon emission process in that the coherent region is extended beyond the 
wavelength. 

The macro-coherent amplification works for both PSR and RENP. Its principle may be 
stated in simple terms as follows [15]. Coherent emission of particles (photons and neutri- 
nos) from a collective body of target atoms is characterized by quantum mechanical rate 
(probability per unit time, not necessarily time independent) which is a squared quantity of 
the sum of amplitudes from many atoms. When plane wave functions of emitted particles 
are extracted, the rate from the collective body is proportional to 

|^e^E,A:.(r-.-.)^^(^-.^^)|2^ (9) 

L 

where ki are momenta of emitted particles including those of neutrinos and rL is the atomic 
position. The atomic amplitude part ALir,t) is expected to be slowly varying with fi in 
the wavelength scale of fej. Unlike the incoherent decay in which the phases of ^^.(r, t) 
from different atoms are random (resulting in the summed quantity Eq. 9 of order A^) , the 
coherent process requires a high level of phase coherence of this quantity to give order A^^ 
to Eq. 9. When a single photon is involved as the only emitted particle of the process, 
ki = k, the maximal coherence region is limited by the wavelength ~ l//c, the inverse of 
the wave number. Nonetheless the coherent amplification oc N'^ {N is the number of atoms 
within the wavelength limited coherent volume) leads to an explosive collective decay as 
first discussed in the celebrated paper of Dicke and later confirmed experimentally (see the 
next chapter on more of this). A typical effect of Dicke super-radiance is the sudden de- 
excitation of all atoms in the wavelength limited coherent volume, emitted photons confined 
in a narrow axial direction. On the other hand, the wavelength limitation is removed when 
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more than two particles are involved as in PSR and RENP and the momentum conservation 
J2i h = holds: the coherent volume may become truly macroscopic without the wavelength 
limitation [15]. In the case of PSR the momentum and the energy conservation limits emitted 
two photons to back-to-back direction and of equal energy eeg/2, if no other phase memory 
is present. The termination of macro-coherent process may be much more rapidly expedited 
than in the Dicke case. Both super-radiance and macro-coherent amplification is a highly 
dynamical process, and one cannot describe its principal feature by a time constant rate. 
One needs a time evolving dynamical equation, which we call the master equation. 

The master equation that describes PSR events including the spatial grating effects (macro- 
scopic polarization varying with the wavelength) has been derived in [14]. The equation is 
given in terms of the Bloch vector components Ri{x,t) ,i = 1,2,3 for medium polarization 
=b ii?2) and the population difference (-R3), and field mode envelopes Ei{x,t) {i = R,L 
denoting right- and left-moving modes) . The single spatial direction, x-axis, is selected as the 
direction of trigger irradiation. Our master equation contains relaxation effects described by 
three time constants, Ts,T2,Ti where the phase decoherence times, T2 for spatially homo- 
geneous modes and T3 for spatial grating modes (both <^Ti) are more important. The 
population decay time Ti may effectively be taken infinitely large in our problem. It is con- 
venient to rescale these and length/time variables x,t using dimensionless variables. The 
dimensionless Bloch vector is rescaled by dividing the excited target number density n, the 
field envelope strength by eegU, while the length/time by i* = 2 / {eegO-gen) ■ Here age oc the 
product of two transition dipole moments. In the double limits of large relaxation times and 
large target length, explosive PSR events may occur with an effective rate proportional to 
the stored energy in \e) / time duration. Dependence of an effective PSR rate on the num- 
ber density is then automatically oc n? and the macroscopic target volume is the relevant 
coherent volume. Thus, if the target length ^ ct^ and relaxation time S> U, one may expect 
explosive PSR, a new phenomenon which may find interesting applications, for instance in 
quantum information. 

Explosive events are not the only important outcome of PSR phenomenon. It turns out 
that static condensate remains after PSR emission and we expect that these condensates 
are described as steady state solutions of our master equation by taking vanishing time 
derivatives, namely static solutions. If one neglects the spatial grating effect, these states are 
expected to become aggregate of many absolutely stable solitons as given in [14]. The most 
important aspect of these condensate states is their stability against two photon emission, 
and their instability for RENP. By its stability the condensate formation makes the signal 
to the background ratio RENP/PSR large. The condensate formation of large active region 
for RENP is then the ideal target state for RENP experiments. PSR (or rather two photon 
emission) has dual roles of importance to realization of RENP; first as a trigger to expedite 
RENP and second as a background of RENP to be rejected. We shall explain how these 
conflicting aspects are reconciled by condensate formation having a very small leakage flux 
at target ends. 
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1.4. More on RENP and PSR 

Let us clarify further the mechanism of macro-coherent RENP amphfication. The electroweak 
amphtude of RENP for a single atom is [8] 

H-yi/ = — ^ ~" — , = Ug^Ugj — Sij/2 , (10) 

where S and d are spin and dipole operators for atomic electron, and U is the unitary 
matrix relating the neutrino flavor, for instance Ug, to the mass eigenstate fj, containing 
mixing angles and Majorana CP phases a, /3 like Ue2 ^ e*" , ^^e3 ^ e''^^~^\ This effective 
hamiltonian is derived in the second order of perturbation of the electroweak theory where 
both Z and W exchange diagrams of weak interaction are involved [8] . Three atomic levels arc 
involved; |e) for the initial, \p) for the virtual intermediate, and 1^) for the final states. Using 
the terminology of atomic physics, one would say that RENP involves Ml x El transition. 
Assume momentarily (proved by macro-coherence) both the energy and the momentum 
conservation of 3-body RENP process. This determines six threshold energies of neutrino 
pair emission at the photon energy lo = ojij = tegfl — {rrii + m.j)'^ j (2eeg). The threshold rises 
of rates are determined by the elements of the neutrino mass matrix, including ^13 and 
Majorana phases (a, /?), as illustrated in Fig.2. 
Let us explain how the Majorana/Dirac distinction comes out. ^ 
The Majorana field [11], [9] can be decomposed in terms of plane wave modes as 

^"^{x.t) = {u{p)e-'^^'+'P-%ip) + u'^{p)e'^^'-'P^{p)) , (11) 

i,p 

where the annihilation bi{'p) and creation bjip) operators of the same type appears in the 
expansion (the index i gives the i— th neutrino of mass rrii, and the helicity summation is 
suppressed for simplicity). The concrete form of the 2-component conjugate wave function 
■u'^ oc i(T2«* is given in [9]. The Dirac case is different involving different type of operators 
bi{p) and dl{p): 

^IjD^^^ t) = Y^ (u{p)e-'^'^+'P-%{p) + v{p)e'^'^-^^-^d\{p)^ . (12) 

i,p 

Neutrino pair emission amplitude of modes ipi,jp2 contains two terms in the case of 
Majorana particle: 

bib] iaiju*{pi)u%p2) - ajiu*{p2)u%pi)) , (13) 

and its rate involves 

-| {aijU*{pi)u''{p2) - ajiU*{p2)u''{pi)f 

= ^|a,,f (|V'(1, 2)|2 + \^p{2, - K(a,,)' (^(1, 2)V(2, 1)*) , (14) 

where the relation aji = (flij)* is used and tp{l, 2) = u*{pi)u'^{p2)- Result of the helicity sum 
^ {tp{l, 2)V'(2, 1)*) is given in [9], which then gives the interference term oc 3f?(aij)^. The first 
term oc la^jp is common to the Dirac and the Majorana neutrino. 



^ In what follows we explain how the interference term arises for the pair emission of Majorana 
fermions using the two-component formalism of [9]. In Sec. 3 a similar derivation of the interference 
term is given using the four-component formalism of the constrained, self-conjugate field = tp. 
Two methods give identical results. 
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The macro-coherence is developed by irradiation of two trigger lasers of frequencies loi with 
the relation uji + (jj2 = ^eg- The development of macro-coherence works for any combination 
{u!i,0J2) and not restricted to oji = UJ2- The frequency at the red side uj\ (< eeg/2) is set for 
detected RENP photon energy. The effective interaction of two fields Ei of frequency ooi with 
atoms is given by 

_ {d)gp-Ei{m)pe-E2 _ ^ ^ 
H27 = = {Ei)i{age)ij{E2)j , (15j 

where m = geS/2me is the magnetic dipole operator, S being the electron spin operator. 
ictge)ij is a tcnsor giving the interaction strength of PSR [14] [16]. The added hamiltonian 
Ti2'y + y.'yu {E = El and w = wi in Eq. 10 ) describes the atom-field-neutrino interaction. We 
shall treat effects of ^27 non-perturbatively, solving this part exactly by numerical means, 
and T-L'yi, is treated in the first non-trivial order of perturbation. 

Two trigger laser irradiation is designed for efficient coherence development. Depending 
on the magnitude of the product of dipolcs that appears in Eq. 15, there may or may not be 
significant PSR emission. What is important for RENP is the later stage after PSR related 
activities. The asymptotic state of fields and target atoms in the latest stage of trigger 
irradiation is described by static solutions of the master equation for time evolution. In 
many cases there is a remnant state consisting of field condensates accompanied with a large 
coherent medium polarization. In the limit of small decoherence time (relaxation time 
for grating modes) this condensate is expected to be identical to a soliton discovered in [14] 
or their aggregate. In any event the asymptotic target state is stable against two-photon 
emission, but RENP occurs from any point in the target. 

Laser irradiation is continued until ~ several times the relaxation time T2 (non-grating spa- 
tially homogeneous modes) of phase coherence and terminated there. This cycle is repeated 
to accumulate detectable level of RENP photons. During a cycle RENP photon is emitted 
within the whole space region within the target length because of the instability of conden- 
sates against RENP. On the other hand, PSR photons are emitted at two target ends due to 
a leakage flux, and not from the inside of medium, because condensates inside the target do 
not emit PSR photons. This way the signal to the background ratio RENP/PSR is largely 
enhanced by the small leakage energy flux at two ends due to QED two-photon process. 

The important dynamical factor r]u{t) in Eq. 8 is given by the bulk integral over all 
target atoms/molecules of the quantity, the absolute magnitude squared of macroscopic 
polarization x the total field strength, both in the dimensionlcss units, as fully explained in 
Sec. 3. The time dependence of r]u{t) disappears in the asymptotic time limit, and this limit 
is described by the static solution of our master equation. Non-trivial static limit exists and 
the asymptotic state consists of the field condensate supported by a macroscopic medium 
polarization. The state is identified as the soliton of [14] in some parameter limit. 

Besides the large bulk/edge rate enhancement there exists a method to selectively detect 
RENP against two photon emission (namely weak perturbative PSR). The magnetic field 
may be used to verify parity violating (PV) effects such as correlated emission of photons to 
the field axis, effect intrinsic to the weak process. Evidently PV effects such as emergence of 
circular polarization and angular correlation of emitted photons to atomic spin are critical 
to prove that the process involves weak interaction. 
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We mention an important point for the target choice of RENP. There is a technical, but 
important reason why one has to look for heavy atoms or heavy molecules as candidate 
targets. This is explained for J = 2 ^ (whose example is the Xe case) and J = — > 
(Yb case) transitions in the following. The Mix El atomic transitions for RENP go through 
the angular momentum change of J = 2 — t- 1 and J=l— >0, orJ = 0— t-I and J = 1 — > 0. 
The Ml transition is governed by atomic matrix elements of (J± ll^lJ). It can be shown 
that these matrix elements vanish in the limit of exact LS coupling scheme. The LS cou- 
pling scheme, however, breaks down for heavy atoms/molecules [17], and the overtaking jj 
coupling scheme evades this constraint, which explains a possible large Ml matrix elements. 

1.5. Relation to cosmology and outlook 

Oni RENP measurements arc sensitive to Majorana CP phases in the following combinations 
in the photon energy thresholds [10]; 



(12) 


(13) 


(23) 


^12*12^13 COs2q; 


C12C13S13 COS 2(^-5) 


S12C13S13 cos2{a - P + S) 



It would be of considerable interest to compare this with the following combination of the 
lepton asymmetry that often appears in lepto-genesis theory [6] : 



, , -2( — )'^sfosin2(a-/3 + (5) + — sin(2a) . (16) 
4n V rn2 '^'^ ^ ^ m2 J 

While the overall factor yi, the Yukawa coupling related to the heavy Majorana fermion, 
is unknown, the other combination of neutrino mass parameters has strong correlation to 
experimentally measurable quantities in RENP. 

If one succeeds in measurements of RENP, one may hope and proceed to detect relic neu- 
trinos of 1.9 K by using the spectrum distortion due to the Pauli blocking effect caused by 
the relic sea of cosmic neutrinos [18]. 



Our experimental strategy towards precision neutrino mass spectroscopy is first to prove 
the macro-coherence principle by QED process of stronger El x El, namely the discovery of 
explosive PSR, and then to control PSR and create ideal form of condensates for preparation 
of RENP. A good target of El x El type is pH2 vibrational transition Xv = 1^0 [14], 
whose molecular properties shall be given in the following sections. 

In the rest of this article we shall give detailed account of theories behind the neutrino 
mass spectroscopy and our current experimental status towards this project. We shall present 
material in a self-contained way even at the risk of overlap with our several past publications 
on this subject. Moreover, the material is presented in order to clarify how we develop 
experimental steps towards the goal of precision neutrino mass spectroscopy. This way we 
give detailed account of PSR and condensate formation so crucial to our neutrino mass 
spectroscopy. 

The rest of the article is organized as follows. In the first two sections we develop theories 
of PSR and RENP based on macro-coherent amplification mechanism. Throughout these 
theoretical sections the natural unit of ^ = c = 1 is used. In Sec. 4 experimental aspects of 
PSR and RENP and the status of our experimental project are described. Several Appendices 
give detailed accounts that supplement the main text. 
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2. Theoretical aspects of paired superradiance 

The twin process of two-photon emission |e) — )• \g) + 7 + 7 is important to experimental real- 
ization of RENP |e) — I^) + 7 + I'ii^j from two reasons; (1) QED proof of the macro-coherent 
amplification mechanism replacing the much weaker neutrino process, and (2) control of this 
QED process to reduce the major background against RENP. In this Section we shall focus on 
theoretical aspects of macro-coherent two-photon emission, termed paired super-radiance, 
or PSR in short. PSR is very interesting and has its own merits in quantum optics. The 
first observation of PSR is our primary goal prior to commencement of RENP experiments. 
Aspects more directly related to real experiments shall be discussed in Sec. 4. 

We deal with atomic processes, either in PSR or in RENP. The fundamental interaction 
of atomic electrons with the transverse electromagnetic field and the neutrino field is briefly 
summarized in Appendix A. The upshot of this approach is to use bound or nearly free 
electron wave functions under the influence of nuclear Coulomb potential for computation 
of atomic transition matrix elements, which prevails in the following discussions. 

It is useful first to compare PSR to the related process of single photon super-radiance [13], 
or SR in short. We shall thus begin with some rudimentary discussions on SR. These two 
phenomena have similarities and differences. Cartoons showing spontaneous emission and 
paired super-radiance are illustrated in Fig. 5 and Fig. 6. At a superficial level this cartoon 
may give an impression that PSR is similar to SR, but the photon energy, the region of 
macro-coherence, and how they are realized is entirely different, as explained below. 

2.1. Super-radiance and extension to two-photon emission process 
2.1.1. Super-radiance. The most common de-excitation process of excited atoms as occurs 
in a dilute gas is the spontaneous decay |e) — )■ \g) + 7, which arises in the presence of a non- 
vanishing electric dipole (El) transition moment between two relevant levels, |e) and 1^). 
The decay follows, as time t increases, the exponential law e~^^ with a decay rate F whose 
inverse is a major portion of lifetime (the inverse lifetime is given by a sum of rates over other 
de-excitation processes). The interaction hamiltonian that appears in the Fermi golden rule 
for the transition rate contains a product of dipole matrix element and electric field E of 
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emitted photon; H = {g\d\e) ■ E. The selection rule for this dipole-allowed transition is the 
parity change and the angular momentum rule A J = ±1,0 (except the strictly forbidden 
J = — > transition) between two states. A typical lifetime would be around 10 ns taking 
atomic energy difference of 1 eV and atomic size of 10~^ cm. Between the same parity 
states the El transition is forbidden and the dominant process of radiative decay may be a 
weaker magnetic dipole (Ml) transition which is caused for instance by an atomic operator 
geS ■ B/{2me) where S is the electron spin operator and B is the magnetic field of emitted 
photon. Ml decay rate is typically a ^ 1/100 smaller than El-allowed decay rate. Besides 
the exponential law the spontaneous decay has an isotropic angular distribution of emitted 
photon unless the initial state is spin polarized (even with spin polarization the anisotropy 
of distribution is a minor effect). 

These features of spontaneous decay are dramatically changed when atoms decay coopera- 
tively. Dicke [13] pointed out in 1954 the possibility of what he termed super-radiance (SR). 
The point here is that for N excited atoms in a volume within the wavelength A of emitted 
photons there is no quantum mechanical way to distinguish a host atom of emitted photon. 
Assuming no decoherence of phases due to interaction etc, all N atoms may decay with 
a definite phase relation preserved from the beginning of the decay process. Suppose that 
several photons emitted in an initial stage have nearly the same direction by a chance, for 
instance along the prolongation axis in a cylindrical configuration of target atoms, then all 
these atoms may undergo cooperative decay. Dicke mentioned a nice analogy of this coherent 
initial state of 1) atoms to an eigenstate of total angular momentum. In this picture 

two relevant states are regarded as spin up and spin down state. Thus, |e) = | t) i Is) = I \)-, 
and the initial inverted state of N atoms is given by an eigenstate of the total angular 
momentum, | J, J) = | t) x I t) x I t) ^ • • • x | |) with J = N/2. The radiative decay of the 
transition \e) ^ \g) is governed by the lowering operator J_ of angular momenta. If any 
phase decoherence does not occur during subsequent radiative decays, then one would arrive 
at a state again given by eigenstate of the total angular momentum, 

\J ,M) = ji-^\J ,J) . (17) 

All these states satisfy the maximal symmetry under interchange of atoms. When M becomes 
0[1] ^ iV, the factor |{J, M — 1| J_|J, M)p, hence the rate of the total de-excitation of the 
collective body, becomes proportional to 0[N'^], as verified by the well known relation, 
J_|J,M) = ^J{J + M){J - M + 1)\J ,M - 1) . The enhancement of decay rate is hence 
due to a/(J + M)(J - M + 1) = 0[iV] for |M| <C A^. The underlying situation is that there 
are 0[N] ways to connect initial and final cigenstates reflecting the indistinguishable nature 
of host atoms. This simple picture breaks down in the presence of de-phasing such as van-der 
Waals interaction, because the maximal symmetry is broken by interaction. 

Two essential features of SR lie in time profile and directionality of emitted photons. 
The temporal feature of the SR dc-excitation as illustrated in Fig. 7 (our own experimental 
result) is the early termination in the form of short pulse at time ~ l/(rA) and duration 
of a short time of the same order. The initial time of cooperative decay, usually called the 
delay time in the literature, may be defined as of this order l/(rA). Distribution of emitted 
direction is uni-directional along the prolongation axis of the cylinder, taken as x direction. 
For simplicity we can ignore transverse field effect parallel to y, z directions in the cylindrical 
configuration. 
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Fig. 7 Super-radiance signals [23] observed in the de-excitation process from Ba 6s6p ^Pi 
to 6s6d ^I?2 state. Note the changes in delay times and peak heights as the density of the 
excited state increases from (a) to (c). 

It would be appropriate now to explain fundamental dynamical variables and their equa- 
tions to further analyze the process. We take the continuum limit of atom distribution within 
the target. The fundamental dynamical variable of medium is then the two-component com- 
plex wave function at each space point ipaix, t),a = e,g, and the field variable is electric field 
E(x,t) (ignoring the vectorial nature). Instead of the wave function it is more convenient to 
use their bilinear forms, called the Bloch vector, R{x,t) = ip*{x,t)aip{x,t) (due to the slow 
atomic motion the norm tp*{x, t)ip{x, t) expressing the atomic number density may be taken 
constant). The Bloch vector is a linear combination of the density matrix element, and it 
obeys in the pure quantum state the von Neumann equation equivalent to the Schroedinger 
equation. What is left to work out is then to write extended von Neumann equation and the 
generalized Maxwell equation incorporating atom-field interaction and effects of relaxation 
in an open system, generally described by two relaxation time constants Tj, z = 1,2 for the 
two-level system [19]. 

Extensive works have been devoted to elucidate the origin of initiation for SR. The com- 
monly accepted view originating from works of [20] is that the inverted state of atoms 
generates quantum fluctuation of medium polarization of 0[y/N], inducing field emission 
due to El coupling in the basic hamiltonian, and it finally develops into a classical coher- 
ence evolution. This gives rise to a natural time scale l/P for initiation. More concretely, the 
origin of quantum fiuctuation lies in the quantum algebraic relation between the medium 
polarization R± and the population difference i?3 that prevails in the entire volume of target: 

[i2+,i?_] =2i23, (18) 

as readily derived from the identity \e){g\g){e\ = |e)(e| etc at each atom. In the completely 
inverted case is of order A'^, and one may take the right hand side of this equation to 
be a c- number of this order. This Heisenberg type of algebraic relation then gives polar- 
ization a quantum fiuctuation of RMS value of = 0[iV], even if its linear average 
value vanishes: {R±) = 0. This 0[^/N] fiuctuation of polarization necessarily couples to field, 
and either zero point fiuctuation of field or spontaneously emitted single photon field may 
stimulate the growth of medium polarization. 
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The classical development at later stages after quantum initiation is well described by a 
non- linear coupled set of partial differential equations [13] for medium polarization R = 
population difference Z, and field envelope E (all functions of time t and x), 



where n is the number density of atoms, and uj is the frequency of emitted photon, while 
d is the dipole moment. This set of equations is often called the Maxwell-Bloch equation. 
The approximations made are (1) rotating wave approximation (RWA) which omits rapidly 
oscillating terms in time, (2) slowly varying envelope approximation (SVEA). The initial 
condition for this set of differential equations is given by random values of initial parameter 
R well described by a Gaussian distribution of width 2-\/N with N the total number of 
participating atoms. The mathematical structure of sine-Gordon equation exists for this set 
of equations, which shall be discussed in Appendix B along with auto-modeling solutions. 

Leaving aside the quantum origin, one can thus clarify the spacetime evolution of SR 
signals, using the classical Maxwell-Bloch equation and taking a Gaussian ensemble of initial 
data set. Solutions give a large fluctuation of the SR delay time, pulse width, pulse height, 
and it is only meaningful to compare SR experiments with theory statistically, and not 
making a shot by shot analysis. Due to the quantum nature of initial stage of SR, the 
phenomenon is sometimes called superfluorescence(SF) instead of SR. 

The super-radiance was first observed in the infrared region of molecular transition [21], 
and subsequently in the optical region, including excited atoms in solids [13]. It was found 
later [22] that the delay time may be shortened by irradiation of trigger along the cylinder 
axis. It is also expected that SR may be controlled to a certain extent. We refer to two 
references in [13] for more comprehensive discussion on SR. 

We show in Fig. 7 our Ba experiment. The global features of numerical solutions along with 
ringing of pulses may be understood by auto-modeling solutions of much simpler ordinary 
differential equation, as explained in Appendix B. 

An interesting and useful point to our RENP project is to view SR as an cooperative 
enhancement of spontaneous decay in the presence of developed coherence. The simple con- 
cept of constant decay rate, namely time independent decay rate per unit time, disappears 
in the phenomenon of super-radiance, because all atoms in the excited state |e) nearly dis- 
appears after the pulse emission and it is meaningless to discuss the latest stage of evolution. 
Another important point to subsequent discussion is that the mechanism of our (triggered) 
PSR is more akin to the triggered SR rather than SF. This makes much easier to handle PSR 
without discussing the stage of quantum initiation as needed for SF. The trigger is essen- 
tial to assist PSR, since the elementary rate of two photon emission involves an effective 
coupling in higher order of QED, and the rate becomes ~ lO^^sec for para-H2Xt; = 1 two- 
photon decay. The quantum fluctuation of medium polarization coupled to field by a much 
weaker effective constant than El is not sufficient for its macroscopic development caused 
by fluctuating fields. This is a fortunate aspect to theoretical analysis of PSR dynamics: 
the semi-classical approach becomes possible in all stages of time evolution for polarization 
development and amplification by trigger. 



dfR = -idEZ , 

dtZ = d{ER* - E*R), 

{dt + dx)E = i2Tru;dnR , 



(19) 



(20) 



(21) 
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Fig. 8 Level structure for two-photon de-excitation. 



2.2. Master equation for paired super-radiance 

We now turn to PSR. Consider three level atom (or molecule) of level energies, > eg > e^, 
as shown in Fig. 8. 

When transition between two lower levels, |e) and \g), is dipole forbidden, the dominant 
de-excitation process of |e) may be two-photon decay |e) — )• \g) + 77, It is described in terms 
of second order, 2x2 matrix hamiltonian [14], 



-Hi 



,E+E- 



aggE+E- 



a. 



Idpedpg 



+ e 



2d^ f 



pe 



^pa 



9,e) 



(22) 
(23) 



(A slightly different notation is used from [14].) Here = ea — Cfe , w = eeg/2 . The upper 
level \p) that has the largest coupling to lower levels, |e) and \g), should be dominant, but 
other competing intermediate states \p) may also contribute to the hamiltonian Hj. In the 
case of trigger for RENP Ml dipole coupling is relevant to one of the two couplings, and 
dpa is replaced by the magnetic dipole moment fipa, typically 0(1/100) smaller. If we ignore 
Lo dependence, the quantities aab,a,,b = e,g coincide with polarizability. For simplicity we 
took isotropic medium and linearly polarized fields, taking E^ as scalar functions E^. The 
diagonal part oc a^a of this hamiltonian describes AC Stark energy shifts, while off-diagonal 
parts oc Oge are effective coupling for two photon emission and absorption. This hamiltonian 
(22) has been derived adopting the Markov approximation by eliminating the amplitude of 
level \p) in the three level {\e),\g), \p)) system. 

The field E is decomposed into positive and negative frequency parts in Eq. 22, and it 
may contain multi-modes, in particular two counter-propagating modes denoted by Er , E^ 
of the same incident frequency, thus E = Er + El. It is useful at this point to recall the 
physical meaning of coupling constants Oab in the interaction hamiltonian. We first note that 
annihilation (a,,) and creation (aj) operators of photon modes are related to complex fields 
by E^ = ai^Juil^V , E^ = a\^Jijj/2V where V is the quantization volume. The important 
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equations obtained after SVEA arc written in terms of envelope functions: 

+age {{Ri - iR2)^^^El + {Ri - ii?2)(+)^^)) , (24) 

+age [{Ri - iR2f^E*ji + {Ri - zi?2)(-^^2)) . (25) 

Quantities R± in the right hand side are medium polarization and i?3 is the population 
difference. The right hand sides of these equations give effects, all in bulk medium, of 
forward scattering oc """^""" n + ""~"gg R^\ backward scattering oc ^""-^ R^^\ RL(right 
left)- pair annihilation oc ageiRi - iR2)^^\ and RR(right right)-, LL(left left) -pair annihi- 
lation oc age{Ri — iR2)^^^ ■ (The pair creation amplitudes appear in conjugate equations to 
those above.) We have used quantities R^"^") e^'^^'^^ as defined by decomposition of three major 
spatially varying components, 

= + i^We^*'^^ + i^S^^e-^''^^ , k = oj>0. (26) 

Quantities i^^^^e^^'^^^ are what are called spatial grating in the literature of non-linear 
optics. The backward scattering terms, and RR,- LL-pair annihilation and creation terms are 
important only in the presence of spatial grating of polarization. Neglect of spatial grating 
is thus equivalent to retaining forward scattering and RL-pair processes, and ignoring all 
other terms whose effects represent propagation effects, as noted in [14]. 

When excitation to the metastable state |e) is done by high quality lasers, there may 
exist a spatial grating in the initial Bloch vector components, which can be incorporated 
as the initial condition in our computation. For instance, excitation by R-moving pump 
laser and R-moving coupling laser may imprint phase factors e^p^^-^) and e~"^<=(*~^) onto 
target atoms, and give the initial spatial grating q-i-{^p-^o)x _ ^-le^^x target atom at the 
position X. This case corresponds to the initial grating of g^^*"^^. But even without initial 
spatial grating two-photon absorption of trigger laser by non-grating components R^p may 
generate spatial grating terms, which thus requires the enlarged set of dynamical variables 
including R^{^\ 

The field equation is to be supplemented by the Bloch equation. In two-photon process 



this reads as 

dtRi = iaee-agg)E+E-R2-iage{e'''^'E+E+ -e-'''^'E~E-)R3-^, (27) 

-1-2 

dtR2 = -{aee-agg)E+E-Ri + ageie'''^'E+E+ + e-'''^'E-E-)Rs-^, (28) 

-'2 

dtR3 = age (i(e^^-*E+£;+ - e-''^^^E-E-)Ri - {e'^^^^E+E+ + e''^'^^ E- E-)R2) - ^ , 

(29) 

E^ = E% + Et. (30) 



Explicit decomposition of these equations into different grating according to Eq. 26 is given 
in [14], using SVEA. We have introduced two relaxation terms inversely proportional to 
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their time constants Tj ,i = 1,2. The origin of these relaxation terms is left mitouched, but 
they may be measured by experimental means. They arise from interaction of subsystem, 
in this case two states in |e), with environment reservoir such as a thermal bath. This 
form is the most general relaxation terms consistent with fundamental principles of quantum 
mechanics in the two-level problem [19]. ^ The phase decoherence time T2 is much smaller 
and more important than the decay time Ti, which may be taken infinitely large for our 
practical purpose. 

It is both useful and convenient to introduce dimensionless quantities and write the 
Maxwell-Bloch equation of our system in terms of these dimensionless quantities. The natu- 
ral time and length scale determined by the field equation is the inverse of a^fewn. We shall 
define this as 

U = . (31) 

A somewhat different notation was used in [14]. The natural field scale is then obtained from 
the Bloch equation as 

E^ = ^=eegn. (32) 

This -E^ is nothing but the energy density stored in the atomic system, when the number 
density in |e) is equal to n. 

It is useful to give numerical values for these quantities. Take the example of para-H2 
molecule and its vibrational two-photon transition v = 1 — >■ of the electronically ground 
state. This example gives 

lO^^cm-^ _ lO^^cm-^ . 2 « 

U ~ 1cm 40ps , ~ 60GWmm-2 ^ . (33 

n n lO^^cm 

Our dimensionless time/length and field units are then (^, t) = (x/t*,t/t*) , Cj = Ei/E^, . 
The Bloch vector components are rescaled by the target number density, to give the dimen- 
sionless quantity = Ri/n. The rs value 1 means the complete inversion in the level |e), 
and -1 the target in the ground state 1^) completely. The medium polarization |ri it zr2| = 1 
implies the maximal coherence. 

We now give the master equation when the spatial grating is present, (the more general 
equation for two color trigger is given in Appendix). In the dimensionless unit they are 

drrf^ = 47_(|eR|2 + |eL|2)rf + 8S>(eiieL)rf + Aj^eRelri'^ + 47_eie^r5+) 

(0) 

-2i(ei - (e^,)2)r(+) - 2i{el - {elf)^-^ - ^ , (34) 

T2 



5.rf ^ = 47_e^elr^"^ - 2^(e^ - {eiry^> + A^-{\eR\' + IclI'V,-'' + 8Q{eReL)rf> - 



(0) 

+2(ei + (ej,)^)r(+) + 2(e| + {elfVi-^ - ^ , (35) 



T2 



' See further below on the relaxation term cx I/T3 when spatial grating terms are present. 
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(36) 

a,rf = -8 (5R(e«e^)rf + 9(eRei)rf + 2z{el - ielf)r[-^ + 2z(ei - {e'k)')r[+^ 

(0) , 1 

-2(ei + (e^,)2)rW - 2(e| + {elfyt^ - '-^^ , (37) 

n 

= 2irf )(e| - {elf) - 2rf (e^^ + (elf) - 8 (5R(e,,eL)r(+) + S>(e«eL)ri+)) - ^ , 

(38) 



(5, - d^)eL = 2 (7+ + 7-rr )eL + 27-rr^e« + - ^^-f )e|^ + '-{r[-'> - zr^^L , 



(39) 
(40) 



7± = . (41) 

Here = Ti/U are relaxation times in the dimensionless unit. In our previous work two 
relaxation terms are taken equal; T3 = T2, but relaxation time T3 for spatial grating terms 
may in general differ from the phase de-coherence time T2 for non-grating terms, as in the 
equations above. 
In the limit of infinite relaxation times the conservation law, 

dr ((ri°V + (rf )2 + (rf )2 + 2r[-\[+^ + 24-\i+^ + 2r(-V5+)) = , (42) 

holds at any point in the target. Without spatial grating terms, another important conserva- 
tion law of energy exists in the Ti — >■ 00 (but with finite T2, T^) limit. In terms of quantities 
of physical dimensions, it is 



I £ dx + 2{\En\' + I^LP)) = -2[\En\' - \E,f'^=o ■ (43) 



This expresses that exchange of energy between field and medium is balanced by a net energy 
flux at two ends of target. This conservation no longer holds in the presence of spatial grating 
terms. But correction due to this violation is expected small because the spatially varying 
factor e^^*'^^ averages out in the spatial integral above. 

Finally, we note that the master equation for PSR has the semi-classical nature: PSR phe- 
nomena can occur without the quantum initiation unlike SR if the trigger field is irradiated. 
The master equation may be solved numerically Tinder a variety of initial and boundary 
conditions. The standard initial-boundary condition (IBC) we use in the present work is 
the symmetric trigger irradiation in which continuous wave (CW) trigger lasers of counter- 
propagating directions and the equal frequency are irradiated. One can think of a few 
interesting cases of IBC for the Bloch vector components such as a complete inversion in |e), 
and a large phase coherence between the two levels, |e) and \g). It turns out that these two 
cases give quite different results. 
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Fig. 9 Time-evolving PSR output fluxes at two target ends, resulting from the symmetric 
trigger irradiation of the power 1 mWmm^-^ under the conditions of the target number 
density n = 1 x lO^^cm"'^, target length = 100cm, relaxation times T2 = Ts = 10, Ti = 10^ 
ns's, and the initial polarization, = 1, all other rj = (no initial grating assumed). About 
40 % of the stored energy in the upper level is released from two ends in this example. 
Vibrational transition Xv = 1 — )• Xv = of pH2 is considered. Right- and left-moving fluxes 
at two ends are identical. 



2. 3. Dynamics of PSR 

For numerical simulations in this section, we have in mind two types of targets as exemplified 
by pH2 and Xe, which has El x El and Ml x El. These target atoms are candidates for PSR 
and RENP. Their coupling matrices (oafe) ^ind other important parameters are as follows: 

PH2; ee3 = 0.52eV (44) 
, , / 0.87 0.055 \ ^^_24 3 1 1 

("<J6) = nr... n on ^m^ 



0.055 0.80 J ' U 1cm lO^icm-s ' 

Xe; \e) = 5p'CP^/2)Qs^[3/2]2, \g) = 5p''So, |p) = 5/(2p3/2)6s2[3/2]i , (45) 
€eg = 8.3153eV , Cpe = 0.1212eV , jpg = 0.28GHz , 

, , / -2.6 X 10"^ 0.47 X 10-3 \ __24 , 1 1 n 
^ ^ \^ 0.47 X 10-3 6.9 J ' 8cm 102icm-3 

Coupling parameters aab for pH2 are given by the polarizability which is justified when 
w <C eeg [14]. 

Results for fields and Bloch vector components have been numerically computed. Previous 
computations [14] have been done without taking into account effects of spatial grating. 

We first compare in Fig. 9 and Fig. 10 results of two calculations; one with spatial grating 
modes included as dynamical variables and another without them, in both cases assuming 
no initial grating modes. Results show a quantitative difference, but qualitatively they show 
similar behaviors of time evolution. In any event it seems difficult to experimentally resolve 
this level of temporal features. All the rest of illustrated outputs are results including grating 
modes as dynamical variables. 

We present in Fig. 9 ~ Fig. 13 (see also Fig. 4 in Sec.l), numerical results for pH2, including 

(i) (i) (i) 

the spatial grating effect. The quantity = r\ + zrg is defined for the total polarization 
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Fig. 10 Flux outputs when spatial grating terms are absent, using tiie same parameter 
set of Fig. 9 
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Fig. 11 Time-evolving PSR right moving output (left moving flux negligibly small), result- 
ing from the symmetric trigger irradiation of the power 1 mWmm~^ under the conditions 
of the target number density n = 1 x lO^^cm"^, target length = 100cm, relaxation times 
T2 = T3 = 10, Ti = 10^ ns's, and the initial polarization of fully grating mode, = 1, all 
other ri = 0. About 30 % of the stored energy is released in this case. Vibrational transition 
Xv = 1 ^ Xv = of pH2 is considered. 

of initial state in these computations. Results of different initial conditions for the spatial 
grating are shown; (1) case of full spatial grating without homogeneous component of the 
coherence, (2) case of null spatial grating, (3) mixture of grating and homogeneous com- 
ponents. With full grating and without any homogeneous component the output shows the 
greatest left-right (LR) asymmetry, a one-sided flux despite of LR symmetric trigger irradi- 
ation. The directionality of the output flux satisfies the phase matching or the momentum 
conservation assuming that the imprinted grating effectively gives medium a momentum. 
The excitation laser is expected to imprint a certain level of spatial grating, and it would be 
important to control this grating for determining the asymmetry of outputs. 

The space-time evolution of Bloch vector components are shown in Fig. 13. It is striking 
that even after dephasing times T2, T3 interesting spatial structures remain finite, suggesting 
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Fig. 12 Time-evolving PSR output fluxes at the left target end in the left panel and 
at the right end in the right panel, resulting from the symmetric trigger irradiation at two 
target ends of the power 1 mWmm~^ under the conditions of the target number density 
n = 1 X lO^^cm"^, target length = 100cm, relaxation times T2 = T3 = 10, Ti = 10^ ns's, and 
the initial polarization, ^ = 0.9, ^ =0.1, all other = 0. About 20% from the right and 
10% from the left of the stored energy are released. Vibrational transition Xv = 1 — )• Xv = 
of pH2 is considered. 




Fig. 13 Space-time diagram of grating and non-grating Bloch vector components 
corresponding to the parameter choice of Fig. 11. 



that explosive PSR event alone is not the whole story of macro-coherent two-photon process. 
Below we shall discuss static or steady state solutions of the Maxwell-Bloch system. 

Summarizing outputs, we may say that effects of grating modes are large when initial 
grating modes are present, but when they are absent, differences are minor. 

The local quantity defined by \E~{Ri — «i22)P/(4eegn^) at each site within the target 
expresses a potentiality of PSR emission, which shall be reflected in fluxes at two ends in 
later times. The spatial integration of this local quantity is denoted by rji^it), which gives an 
instantaneous potentiality of PSR activity. This global function is illustrated in Fig. 14 up 
to a time of order T2. The time duration of the plateau region roughly gives the duration of 
explosive PSR emission. 
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t (ns) 

Fig. 14 Potentiality of PSR emission given by the spatial integral of \E~{Ri — 
ii?2)P/(4eepn^) corresponding to the parameter set of Fig.9. 

In Sec. 4 wc show the case of weak PSR outputs in which the output flux is linearly 
proportional to the trigger flux. There is a set of threshold parameters that divides the 
explosive and the linear PSR events. But due to a high level of non-linearity it is difficult to 
locate this threshold exactly. 

The results so far are outputs for pH2. Xe results have different aspects, since their 
two-photon coupling is caused by Mix El, much weaker than ElxEl of pH2 vibrational 
transition. 



3. Theory of macro-coherent radiative emission of neutrino pair (RENP) 

In this section, we present a theory of RENP [8] , [9] and numerical results based on theoretical 
formulas thus derived. The RENP process is treated in the second-order perturbation of QED 
and the weak four Fermi interaction, while the time evolution of the electromagnetic field in 
the target medium is dealt with non-perturbatively in the semi-classical approximation. A 
Feynman-like diagram of the elementary RENP process is depicted in Fig. 15, in which the 
nuclear Coulomb interaction is emphasized. The non-perturbative part requires numerical 
solutions of the master equation of Maxwell-Bloch type under two color trigger irradiation. 

In Appendix A a brief description of electroweak interaction of atomic electron is given, 
which is the basis of our following computations. 

3.1. Coherent neutrino pair emission from atoms / molecules 

We begin with a description of the elementary process |e) \g) + ^{k) + i^iip) + i^j{p'), in 
which an electron bounded in an atom of the macroscopic target undergoes the QED and 
weak interactions. The amplitude is expressed as 

^. ^ iMoM^ e-(^>^^+^^')--(27r)(^(6e, - u - - E^,) , (46) 
epg - u 

where epg = ep — eg, to = \k\ is the energy of the photon, M.d{w) denotes the El (weak) 
matrix element and the possible sum over the intermediate state \p) is to be understood. We 
assume that the transition between |e) and \p) is the Ml type responsible to the neutrino pair 
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Fig. 15 RENP process. The many wavy hnes between e and represent frequent 
Coulomb interaction between atomic electron and nucleus. 

emission and that between \p) and \g) is the El type responsible for the photon emission. 
We explicitly show the dependence on the position of the target atom Xa in Eq. 46. 

For an ensemble of atoms, the total amplitude is the superposition of the contributions from 
all the atoms and proportional to exp [—i{k + p + p') ■ Xa] - For a macroscopic ensemble of 
A'' atoms with a volume V, this summation approximates to {N/V){2Tr)^6^{k + p + p'), which 
means the momentum conservation among the emitted photon and the neutrino pair in addi- 
tion to the energy conservation represented by the 5 function in Eq. 46. Thus, the kinematic 
configuration that satisfies the momentum conservation acquires an extra enhancement in 
the rate. This is the essence of the macro-coherent amplification mechanism. 

Due to the energy-momentum conservation, the RENP process with the macro-coherence 
resembles a relativistic three-body decay of the parent particle of mass €eg into a photon 
and a pair of neutrinos. However, the energy and the momentum of the photon are dictated 
by the classical radiation field induced by the trigger laser. Therefore, the emitted photon 
has no phase space and we do not have to introduce the three-body phase space. Thus, the 
differential rate of the process for a macroscopic ensemble of atoms is written as 

dT, n^V^p^^^d^2, (47) 

where n is the number density of atoms, d^2 represents the Lorentz invariant two-body 
phase space of the neutrino pair, and the summations over the atomic spin states and the 
neutrino helicities are implicit. The explicit form of d^2 is given by 

.#. = (2.)V(,-p-p')p||^p^, (48) 

where Ep(,) = ■\J^'i(^j) + P^'^^ is the neutrino energy of mass "^i(j) and q'^ = {eeg — oj, —k) 
represents the four momentum of the neutrino pair. 

The El matrix element is defined as = —{g\d\p) ■ E with E being the radiation field 
developed in the target by the trigger laser irradiation, while the weak matrix element is 
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expressed by 

M% = %{<P, ^)i^j{P, A')| Yl - l^HmvabJi^ - aabJ'^) , (49) 

a,b 

where the atomic vector (axial-vector) current J(/(J^) is defined by Jy = (5167^61^) (J^ = 
{g\e^^j5e\p)). Coefficients Vab and Uab are given in terms of elements of neutrino mass matrix, 
as given by Eq. A12 of Appendix A. In the non-relativistic approximation for the atomic 
wave functions, ~ and ~ (0, 2{p\S\e)) in the rest frame of the atomic system. Hence, 
we neglect Jy hereafter. Evaluating the neutrino matrix element, we obtain 

M% = -^ (a,,^^ - SMajiE^^ J a, , (50) 

where 



(51) 



and Sm = 0(1) for Dirac(Majorana) neutrinos. The left-handed current Lfj, which comes 
from {a,b) = in Eq. 49, is common for both Dirac and Majorana cases; while the 

right-handed current i?^- from (a, b) = (j, i) with the use of the relation between the four- 
component spinor, Cv?' = v, is allowed only for the latter case. Technically, this chirality 
flip is understood as a property of the charge conjugation operation in the Dirac theory. 

The square of the weak amplitude Eq. 50 consists of four terms: Lf-Lj^, Rf-R,]'^, Lf-R^'^ 
and R^-L]^. The sum over the neutrino hehcities and momenta leads to 



•' A,A' A,A' 



A. 



_Ki 



A, 



^3 



3 1 



+ 2(1 + 



mj + rrij 



(52) 



j d^, ^r,< = / rf^2 E = -^mm.g^'^ , (53) 



A,A' " A,A' 



where 



Al = l-2^^+'—^, (54) 

and = €eg{eeg — ^) is the invariant mass of the neutrino pair. The cross terms Eq. 53 
exhibit the identical particle effect inherent in Majorana neutrinos. They arise from the 
interference of two currents of different chirality and R^- defined in Eq. 51 and disappear 
in the massless limit as it should. 
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Thus, taking the statistical factor for the Majorana case into account, tlie single photon 
spectral rate at frequency uj from the position x of unpolarized targets is given by 



,^ , , GlVdx 



E {uj,x,t) 



Ri{x, t) — iR2{x, t) 



Ceg{u)Fiu;) 



r . ._v^ (g|rf|p)-(y|rf|g)(e|%)-(p|5|e) 



(55) 
(56) 



F(a;) = Aij{BijIij{Lo) - 5MBf^ mimj)d{uij - w) , B^j = [a^.f , = ^{aD , (57) 



2- 



m? + m3 (m? 



m 



2\2 



+ 



mf + (mf — m^Y 

1 + " o ^ - 2 ^ ' 



, (58) 



where E^{uj,x,t) = + -^l negative energy part of right- and left-moving fields 

within the target developed from trigger laser of frequency uj < eeg/2 and is limited in its 
magnitude by the maximum energy stored in the target, \E\'^ < eegU. 
We note that the rate dTj2u is proportional to the space-time dependent factor, 

2 

(59) 



iRi 



which signifies how much the field-medium system has been activated in favor of RENP. 
It is thus important to obtain this factor as large as possible in the target. We define the 
following dimensionless quantity to express the field-medium activity for RENP, 



»7a;(4,r) 



1 



E' 



Ri - iR2 



ri — ir2 



(60) 



where (C,t) = {x/t*,t/t*) with given in Eq.31, the slowly varying dimensionless envelope 
functions eR^L{i,T) are introduced as 



(61) 



and the dimensionless Bloch vector components ri(^,r) are defined by Ri = nvi. Further, 
using the grating modes = r-°^ -|- r^+^e^*'^**^ -|- r^~^ e~^'^^^*^ , we obtain 



(62) 



where 



„(+)|2 



»|2 



(0)* (+) ^ JO) J-)* 



with '^^ = r^°'^^ -|- ir2°'^^ . We expect that the 77 factor weakly depends on the frequency 
a; as is described below. 

The formula here gives rate of a kind of stimulated emission from the developed trigger 
laser of field strength \E\'^. Although both the energy and the momentum is conserved at 
the level of elementary process, the phase space factor is different from a typical electron 
spectral shape of three body decay into almost massless particles such as the muon decay, 
since the quantity \E\'^ here is replaced by a kinematic factor oc co (electron energy) in the 
muon decay. Thus, the photon spectral rate as a function of uj is markedly different in 
the small energy region from that of muon decay. Near the neutrino pair threshold region 
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the spectrum resembles that of muon decay, although the difference reflects effects of finite 
neutrino masses. 

As stated in Section 1, RENP may occur at any point of target since even condensate or 
soliton formation, capable of blocking PSR, cannot stop RENP. RENP is a bulk process, 
while PSR occurs only at target ends after soliton formation. The total RENP rate observed 
at both ends of the target is given by integrating Eq. (55) over x. As a result, assuming the 
sum over possible intermediate states is dominated by a single state \p) , we obtain 

T^2Au;,t) = ToI{u)i]^{t), (64) 
To f,'^^''^\ 2 J, + 1)C.,, (65) 

where I{io) = F{ui)/{epg — w)^, and denotes the El transition rate from \p) to \g). The 
atomic spin factor Cgp is defined by 

Y,{pMp\S\eMe) ■ {eMe\S\pM;) = Sn^M^Cep , (66) 

^ Me 

and Cep = 2/3 for Xe atom. The dynamical factor r)^{t) is given by 

vUt) = vSit) + r,t{t), (67) 

where 

f fL/t, 

"^^^^^ = 41/0 " + ^^^^ ' ^^^^ 



and 



f rL/U 



(69) 



correspond to the rates at the right and left ends respectively. The overall rate Vq has a 
dimension in 1/time, while and ?/i^(t) are dimension-less. 

A nice feature of the formula Eq. 64 is that one can separate the spectral feature given 
by I{oj) that serves to determine neutrino properties from the absolute rate Fq and the 
dynamical factor ri^{t). The time dependent function rjujit) is sensitive to how the initial 
coherence between |e) and \g) has been prepared and to the quality of trigger laser. Indeed, 
calculation of 77^, (t) requires numerical solutions of the master equation. We thus discuss two 
issues, the spectral shape and the overall dynamical rate factor, separately in the following. 

The limit of massless neutrinos gives the spectral form, 

/(c.;m. = 0)= i2(e,,-^)^ ' ^"'^ 

where the prefactor of Bij = 3/4 is calculated using the unitarity of the neutrino mixing 
matrix. On the other hand, near the threshold these functions has the behavior of Ajj(a;) 

OC -^UJij — OJ. 

In what follows we first calculate numerically (analytic results are of little use to experi- 
mental design) the spectral shape for two target candidates, Xe atom and I2 molecule, 
which greatly differ in their available energies, 8.32 eV (Xe) and 0.88 eV (I2 relevant vibra- 
tional transition). These two targets have been selected for their large Mix El couplings and 
their large difference in closeness to the minimum expected mass of neutrino heaviest pair 
~ O.leV. This makes it possible to determine the relative easiness of undetermined neutrino 
parameters in realistic experiments. 
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Table 1 The threshold weight Bij = lajjf = \U*^ej - Sij/2\'^. 





B22 


^33 


B12 + B21 


B23 + -B32 


B31 + -B13 


{42^3 - 1/2)2 


{sj.cl, - 1/2)2 


(sh - 1/2)2 


ZC]^2*12'-13 


o„2 2 2 
^*12'-13'^13 


0^2 2 2 

'^''12''13*13 


0.0311 


0.0401 


0.227 


0.405 


0.0144 


0.0325 



5.^. Neutrino properties extractable from the photon spectrum 

Recent reactor and T2K experiments have given compclUng evidence for a sizable ^13, which 
determines relative weights of neutrino pair emission thresholds in our RENP process. Here 
we summarize the neutrino parameters employed in the following numerical calculation: ^ 

Amii = 7.5 X 10"^ eV2 , \Aml^^2)\ = 2-32 x 10"^ eV^ , (71) 

sin^ ^12 = 0.31 , sin^ 013 = 0.025 , sin^ ^23 = 0.42 , (72) 

for the NH(IH) case. It is useful to recall weight factors Table(l) given in reference [12]. These 
absolute weights are applied to the Dirac neutrino pair production, while the rate of the 
Majorana neutrino pair production at the crossed channel ujij, i ^ j contains the additional 
term proportional to B^ in Eq. 57, which depends on Majorana phase combinations, 

cos2q;, cos2(^-^), cos2(a-/3 + J), (73) 

at (1,2), (1,3) and (2,3) thresholds. 

On the other hand, the threshold locations are separated by energy differences eeg/2 — 
LOij = (mj + mj)2/2eec;- Hence the separation becomes larger as the atomic energy €gg is 
smaller for fixed neutrino masses. A numerical illustration of the threshold locations is given 
below for Xe atom in Fig. 17. 

3.2.1. Absolute masses and mass hierarchy. Experiments are easier when detectable rates 
are larger. The single photon spectral rates for macro-coherent RENP are larger for larger 
available energy if dipole and spin matrix elements are the same. Thus, we first calculate 
RENP rate for Xe which has a large energy difference between the first excited state, as seen 
in Fig. 16, which is metastable, and the ground state. The overall rate for Xe is estimated 
as Fo ~ lHz(n/1022cm-3)3(y/102cm3). 

In Fig. 17, we show the threshold locations, that is, the threshold displacements from the 
maximal photon energy for the massless neutrinos, as functions of the lightest neutrino mass. 
The multiplet structure is clear: a triplet, a doublet and a singlet for the NH case and the 
ordering of the multiplets is inverted for the IH case. Since the energy resolution of trigger 
lasers is not worse than O(10~^) eV, the identification of thresholds is possible provided that 
we have a sufficient statistics of data. The expected number of events at each threshold is 
proportional to the corresponding weight factor shown in Table 1. The threshold of (1,2) has 
the largest weight and (3,3) is the second one. Observing these thresholds, we will be able to 
determine the corresponding neutrino masses and examine their consistency with the known 
mass squared differences. 



* These values are slightly different from those quoted in Eq.2. The main conclusions in this section 
should not be affected by these differences. 
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Fig. 16 Xe energy level diagram for RENP experiment. 




mi (meV) (meV) 

Fig. 17 The locations of the six thresholds, the distances in energy from the half of level 
difference €eg/2 given in meV unit, as functions of the smallest neutrino mass for the case of 
Xe atom. The left(right) panel shows the NH(IH) case. The pairs of numbers in parentheses 
indicate the thresholds and two prominent thresholds are printed in the bold font: (1,1) in 
black, (1,2) in blue, (2,2) in green, (1,3) in brown, (2,3) in yellow and (3,3) in red for the 
NH case; (3,3) in black, (1,3) in blue, (2,3) in green, (1,1) in brown, (1,2) in yellow and (2,2) 
in red for the IH case. 



The global spectral feature of RENP from Xe J = 2 excited state 5p5(2p3/2)6s ^[3/2] j=2 
is shown in Fig. 18, while the enlarged threshold region is shown in Fig. 19. The locations 
of the thresholds corresponding to the three values of the smallest neutrino mass mo differ 
substantially. This feature can be used to determine the absolute neutrino mass. 

There are two dominant thresholds and they are supposed to be identified in a relatively 
earlier stage of experiment. The most significant threshold comes from the emission of the 
neutrino pair (1,2) and the second one is due to the (3,3) pair in both the NH and IH cases. 
The relative strength of these two thresholds is a powerful tool to determine the hierarchy 
pattern: The threshold at higher photon energy is stronger than the one at lower energy in 
the NH case. While the former is weaker than the latter in the IH case. 

3.2.2. Dirac or Majorana. Although Xe gives a good chance of determining the overall 
neutrino mass scale or the largest neutrino mass, along with distinction of IH and NH, it 
would be too demanding to differentiate Majorana neutrinos from Dirac neutrinos in Xe 
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Xe NH and IH,mO=20meV 




Fig. 18 Xe RENP spectrum I{oo). In this global plot, the NH and IH spectra completely 
overlap each other, and it is difficult to extract neutrino parameters. 




Fig. 19 Xe RENP spectrum /(w) in the the threshold regions, using values of the smallest 
neutrino mass, 2 and 20 meV, for two cases of IH (in dashed curves) and NH (in solid). 



experiment, since their rate difference due to the interference term is too small in this case. 
We may characterize Xe RENP measurement as the first generation RENP experiment, since 
it serves to identify the process itself along with determination of an important neutrino 
parameter, the largest neutrino mass, although there are better targets for determination of 
other neutrino parameters. 

In a recent paper [12] it was suggested that a smaller available energy of a fraction of eV, 
but still much larger 0.1 eV, might work for Dirac vs Majorana distinction. A candidate target 
is I2 molecule, which has about 1/10 of Xe available energy. The I2 RENP experiment can 
be characterized by the next generation experiment of precision neutrino mass spectroscopy. 

Before presenting RENP spectral shape from I2 molecule, it might be instructive to com- 
ment on molecular energy levels and their associated wave functions [17]. As is well known, 
di-atomic molecules have vibrational and rotational energies associated with relative atomic 
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Table 2 I2 parameters. 



Electronic state of I2 


£'e(eV) 


re(eV) 


re (10~^ cm) 


LOe (cm~^) 




1.556 





2.666 


214.5 


A' 2u{^Uu) 


0.311 


1.245 


3.079 


108.8 


A lu(3n„) 


0.203 


1.353 


3.129 


88.3 



motion besides electronic motions. Electronic, vibrational, and rotational energies are well 
separated due to three well separated velocities. This makes it possible to compute electronic 
energy levels assuming fixed positions of two nuclei, leading to the concept of the potential 
energy curve of each electronic configuration as a function of nuclear distance. Near the equi- 
librium position of nuclei the potential curve may be approximated by a harmonic function 
whose frequency oj^ gives energy levels of harmonic oscillator ujy{v + 1/2). We must deal 
with two potential curves corresponding to two different electronic configurations, both for 
Ml and El transitions. This gives two overlapping integrals of vibrational wave functions. 
Each of these overlaps is called the Franck-Condon (FC) factor and calculated assuming 
relevant vibrational quantum numbers of respective electronic configuration. Some rudimen- 
tary example of these is illustrated in Appendix, using the Morse potential for vibrational 
motion. A work is in progress for accurate computation of relevant potential curves for I2, 
in collaboration with theoretical chemists [24]. 

We have found that three relevant electronic states for RENP are X (electronically ground 
state) for \g), A' (metastable state) for |e), and the largest overlap comes from A {\p)). Their 
experimental data are taken from experimental values tabulated in Table 2. They are also 
used to construct the Morse potential in Appendix. 

Here the energy unit lcm~^ = 1/8065. 73eV is used for the curvature of the potential curve. 
The dissociation energy De is the energy difference between the potential minimum at Tg 
and level at dissociated atoms. Conceptual diagram of relevant potential curves that fit to 
these experimental data is shown in Fig. 20. 

We show a threshold feature of RENP photon spectra for the transition between A'v=l 
and Xv=15 of I2 in Fig. 3. Although the difference in the spectral rates between Dirac and 
Majorana is not so significant even in I2, it is possible to discriminate between them in 
principle. 

3.2.3. CP violating phases. The rates of emission of vi^s and of 1/2^3 are suppressed by 
the weights Bis + -B31 = 0.033 and B23 + B32 = 0.014, respectively. Thus, it is most advan- 
tageous to obtain information on the Majorana phase a at the (1,2) threshold, provided that 
the corresponding interference term is not suppressed by the smallness of the factor mim2. 

We note that a linear combination of the CP violating phases given in Eq. 73 enter into 
the expression for the effective Majorana mass measurable in the neutrinoless double beta 
decay: 

I ^ mii/^ |2 = 5.8 X 10-^ ml + 9.2 x lO^^ ml + 4.5 x 10"^ ml + 4.1 x 10"^ mim2 cos(2a) 

i 

-F3.2 X 10~^ mima cos 2(/3 - 5) + 1.5 x 10~^ m2m3 cos 2{a- 13 + 5), (74) 
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Fig. 20 Conceptual diagram of potential curves of relevant electronic states of I2 for 
RENP; |e) = A' ,\g) = X ,\p) = A. Vibrational levels are suppressed for simplicity. The 
curves are drawn based on [25]. 

where the values given in Eq. 1 are used for the mixing angle parameters. The RENP 
experiment has an advantage of independently determining all three CP factors Eq. 73 
instead of their single linear combination. 

Sensitivity to the CP violating phases in I2 RENP spectrum is not as great as one might 
have expected, but differences for different parameter sets {a, (3) reaches several to 10% for 
the inverted hierarchy case, as seen in Fig. 3 of Section 1. For measurements of CP phases in 
the NH case one should perhaps look for better candidate targets of smaller available energy 
than I2. 

We show here I2 RENP spectrum results for a smaller uiq = 5meV, which gives smaller 
differences, as given in Fig.21. 

3. 3. Estimation of the dynamical RENP factor based on asymptotic solution 
Our numerical outputs on RENP have so far been limited to the spectral shapes, which are 
crucial to determine the neutrino parameters. Some comments on the a; dependence of rji^^t) 
in the fundamental formula Ty2ui^^, t) = ^oIi^)ili^it) are now in order. Numerical simulations 
are necessary to compute this quantity at each photon energy uj, solving the master equation 
Eq. Bll ~ B18 in Appendix B for two-color trigger irradiation. The u> dependence of r]i^{t) 
thus originates from u dependence of basic couplings Oabi^jj) that appears in the master 
equation, in particular through 7_ = {uee — 0!gg)/2age and age- As seen in the formulas for 
these in Appendix B, both of these are smooth functions of oj. Furthermore, if the upper 
level \p) lies much above eeg (the energy difference between |e) and \g)), these aafe('^)'s are 
nearly uj independent. In any event, one expects that variation of r/aj(i) in the vicinity of 
each threshold ujij can be neglected, and the theoretical formula of I{uj) reflects the spectral 
shape of experimental data in the threshold regions. This way one can extract important 
information of neutrino parameters directly from measured data. It is however desirable to 
measure the quantity rj^{t) or their asymptotic value at times ^> ?2, T3 for the overall fitting 
to the entire spectral shape. This method is described in Sec. 4. 
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12 A'v=l->Xv=15: mO=5meV 
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Fig. 21 I2 RENP spectrum between A' v=l and Xv=15. The Majorana (in colored solid 
for NH and in colored dashed for IH) vs Dirac (in black solid for NH and in black dashed 
for IH) cases are compared. A Majorana CP phase combination (7r/2,0) (in red) is taken, 
with the smallest neutrino mass 5 meV. The vertical scale is in arbitrary units. 



Whether one can actually perform RENP experiments with meaningful results depends 
on the absolute magnitude of rate Pq: if it is too small, the spectral shape com- 
putation would be useless. For instance, although the Xe RENP rate Pq is of order 
1 Hz (n/10^'^cm~^)^ (y/lO^cm^), the dynamical factor r]i^{t) might give a further large 
suppression. 

Due to the complexity of our non-linear system our simulations so far are limited to the 
explosive cases and the linear regime, mainly for pH2. We feel that the late time evolution in 
the intermediate output regime (case of large Ml xEl coupling) is very important for RENP, 
which has not been explored extensively due to technical complexity. Moreover, the RENP 
target choice is very limited so far: Xe and I2 (whose absolute rate is under study[24]). The 
case of Xe transition involves a smaller Oge (due to Mix El nature of RENP), but a larger 
7_ ~ —7300 (the corresponding number for pH2 ~ 0.6), which makes numerical simulations 
for r]i_j(t) harder. Under this circumstance our conclusion on the absolute rate is premature 
and we shall further update the simulation. Nonetheless, it would be useful to summarize 
what we have done so far. 

We expect that the final stage of field and target state within the medium is described 
by static solutions of our master equation, which gives a long time behavior and is stable 
beyond the largest relaxation time Ti. These static states are solutions obtained by taking 
vanishing time derivatives of the master equation, Eq. Bll ~ B18 in Appendix B. We first 
solve Bloch vector components, both grating and spatially homogeneous modes, in terms 
fields. These are inserted into the field equations to search for static solutions by solving the 
coupled system of ordinary differential equations. 

Results of spatial profiles of the static asymptotic solutions are illustrated in Fig. 22 and 
Fig. 23 in dimensionless units. The master equation for the single color mode, Eq. 34 ~ 40, 
has been used for simplicity. Due to the choice of a small this solution resembles the soliton 
solution given in [14], but it gives a relatively large tji^ (defined by Eq. 67 ^ 69), too. This 
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Fig. 23 Field intensity and the dynamical factor ij^ for the same parameter as in Fig. 22 
The red (blue) curve shows (e^), and the black line is 7]^^ x 100. 



asymptotic case gives a time independent dynamical factor; ryaj(i) = (constant). Whether 
a large t]^^ factor is realized or not in actual experimental situation needs further study, 
both from theoretical and experimental sides. Moreover, it should be clarified how the time- 
evolving dynamical factor r]^{t) approaches the asymptotic constant rj^^. Also, relevance of the 
field condensate making up the asymptotic rj^j value to solitons in [14] should be investigated. 
(The soliton is distinguished from a mere condensate by the presence of topological quantum 
number.) 
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4. Experimental aspects of PSR and RENP 

4.I. Overview and strategy towards PSR/RENP experiments 

4.1.1. Experimental principle of the neutrino spectroscopy with atoms. We begin this 
Section by reviewing the experimental principle of our neutrino spectroscopy. First of all, we 
prepare a collection of excited state of atoms or molecules. We then observe photon spectrum 
in a radiative emission of neutrino pair (RENP) process denoted by 

\e) ^ \g) +7 + + (75) 

where fj's are neutrino mass eigenstates. Rich physics information exists in the photon 
spectrum, especially near the threshold regions given by 

_ eeg {rrii + mj)'^ 



2eeg 



(76) 



where egg is the energy difference between |e ) and \g), and denotes the neutrino 

masses. The physics observables include 

o neutrino absolute mass scale, 

o mass hierarchy pattern, normal or inverted, 

o mass type, Dirac or Majorana, 

o CP violating phases, a and /3 — S, 

among others. Detailed discussion of the physics objectives can be found in Sec. 3. 

A key notion in our experimental principle is the macro- coherent amplification mechanism. 
In a word, the amplification is based on their collective de-excitation among all of the 
coherent atoms. It enhances not only the rate of the de-excitation process but also guarantees 
a "momentum conservation law" among the emitted particles (photon and two neutrinos in 
the case of RENP). The physics information can be extracted from the photon spectrum 
with the help of this momentum conservation law. Since we crucially rely upon this macro- 
coherent amplification mechanism, we plan to prove its validity by a separate experiment. 
We discuss this proof-of-principle experiment in more detail below. The readers are referred 
to Sec. 2 for detailed discussion on the theory of the macro-coherent amplification. 

Another important feature in our experiment is use of the trigger laser. It initiates coherent 
de-excitation process, and in effect selects the portion of the photon spectrum we like to 
investigate. In this sense, our neutrino spectroscopy uses a "narrow-band spectrometer" . We 
note that energy resolution of the spectrometer is determined by the laser line width, which 
can be much less than a few hundred MHz, or ~ 1 /xeV. Detector energy resolution is not 
critically needed for spectroscopic information, but it helps to reject background signals. The 
role of the trigger laser goes beyond the initiation of coherent de-excitation process; it helps 
developing condensate formation, another key notion in the RENP neutrino spectroscopy. 
The condensate, a system of coherent atoms and fields, is an ideal target state for RENP 
process. 

As discussed in Sec. 3 in detail, the final RENP rate r^2i/ is given by 

r^2A^,t) = ToI{Lj)riUt) (77) 

where Fq is the overall constant rate and I{cj) is the spectrum function containing physics 
information. The third factor rjoj, termed as the dynamical RENP factor, is proportional to 
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the squared product of the field E~ and the medium's coherence — iR2)/2, averaged over 
along the target. In order to extract I{oo), we need to divide experimentally observed rate 
^'y2i/ by To and rj^. The constant Fq can be calculated accurately while the RENP factor rj^ 
may be obtained by detailed simulations. It is, however, highly desirable to measure rj^^ more 
directly. Wc plan to determine it experimentally by measuring the coherence and population 
difference (rs) of the target medium as well as the field intensities at both ends. With aid of 
the master equations, 77^; can be determined from those quantities determined above. One 
method of measuring coherence (ri,r2) for the case of Xe atoms is described in Sec. 4.3. 

We detect single photon emission by an appropriate photon detector. The detector may be 
UV to IR photo sensors, depending upon the value of initial energy gap between metastable 
and ground states, but most likely we must use IR sensors to investigate the detailed nature 
of neutrinos such as Majorana-Dirac distinction and CP phases. Wc expect the photon 
emission direction to be along the trigger laser direction because macro-coherence mechanism 
preferentially develops the process with this mode. Thus detector solid angle coverage need 
not be large. 

4-1-2. Paired Super-Radiance. Exploitation of the macro-coherent amplification mecha- 
nism is an essential feature in our neutrino spectroscopy. Since this is a new phenomenon not 
discussed in the past literature, we should first prove and study it in detail experimentally. 
We like to do it with a QED process called paired super-radiance (PSR), a process which 
emits two photons in the decay from an excited state to the ground state (|e) — >■ l^r) -|- 77). 
If the macro-coherent amplification is in operation, a collection of atoms emit coherently 
a pair of photons with the same energies. In particular, direction of the two photons is 
opposite (back-to-back) when the medium's initial macroscopic polarization (coherence) is 
spatially-homogeneous. Thus observation of such radiations in a similar fashion to SR is 
an unambiguous proof of the principle. There are several advantages of proving the new 
principle with PSR. The most important is its rate: it is expected to be much larger than 
RENP since it is a QED interaction. Characteristic features of the radiation, back-to-back 
with equal energies, make experimental distinction easy from other possible backgrounds. 
The conditions and natures of PSR, mainly explosive PSR, are already shown in Sec. 2. We 
summarize them here briefly, adding some more studies concerning to non-explosive events. 
As shown in Sec. 2, explosive PSR events have the following characteristics: 

o Major part of energy stored initially in medium is released as two photons in an 
extremely short period of time compared with medium's natural lifetime. 

o There exists certain threshold conditions above which explosive events take place. 

o The threshold conditions are governed predominantly by target density and length, 
initial coherence and de-phasing time, etc.. 

o Left-right symmetric (back-to-back) events occur when initial coherence is spatially- 
homogeneous (i.e. r^i '' or r^^ ~ 1 ) and trigger laser is injected from both ends. 

o Asymmetry between left and right flux grows as a larger initial value is assumed for the 
spatially-inhomogeneous coherence. 

It is of great importance, from both theoretical and experimental view points, to study what 
is expected if some of the threshold conditions are not fulfilled. Below we show some results 
obtained by numerical studies performed with these motivations. The process of interest is 
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Fig. 24 Right-moving output flux from pH2 vibrational excitation level {Xv = 1). The 
initial coherence is assumed to be spatially homogeneous (r[^^ = !)• The black, blue and red 
lines correspond to trigger laser intensities of 1 mW/mm-^, 10 mW/mm-^, and 50 mW/mm^. 
Note logarithmic scale for the vertical axis. Left flux is identical to the right flux. 



the same process treated in Sec. 2; de-excitation process of para hydrogen molecules from a 
vibrationally-excited state {Xv = 1) to the ground state {Xv = 0). Figure 24 shows right- 
moving output flux. The three different curves in black, blue and red, correspond to three 
different input trigger laser intensities of 1 mW/mm^, 10 mW/mm^, and 50 mW/mm^. 
Other important parameters are 

n = 8 X 10^° [cm~^], L = 10 [cm], rf^ = 1, Ta = = 0.2 [ns] (78) 

where rf^ means initial coherence of spatial homogeneous mode, with all other coherence set 
to zero. The important fact to these results is the values of T2 and T3; here a factor 50 smaller 
value is assumed than the typical values employed in Sec. 2. As shown in Fig. 24, events are 
found to have much simpler time structure compared with explosive events shown in Sec. 
2. Also their peak heights are much smaller and proportional to the input trigger intensity. 
Although these events appear less dramatic, if observed, they bear great significance because 
the two-photon process, which has the natural life time of ~ 10^^ sec, is seen in a time period 
of nsec. 

Initial coherence of spatial grating modes, r^2; affects substantially partition between left 
and right output flux, as already pointed out in Sec. 2. Figure 25 shows an output flux when 
the initial coherence is set to 

r[°) = 0.1, ri+) = 0.9, (79) 

with all other parameters left same as in Fig. 24. In this case, the right output flux is much 
smaller than the left. The fact remains true that the output flux is proportional to the input 
laser intensity. 

The requirements for the PSR target, different from those for RENP in certain aspects, 
can be summarized as follows. 
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Fig. 25 Left and right moving output flux from pH2 vibrational excitation level {Xv = 1). 
The initial coherence assumed is r^'' =0.1 and r^^^ =0.9. The black, blue and red lines 
correspond to trigger laser intensities of 1 mW/mm^, 10 mW/mm^, and 50 mW/mm^. Note 
logarithmic scale for the vertical axis. 



o The initial state must be metastable (El forbidden), and must have at least one inter- 
mediate state of higher energy connecting between initial and final (ground) states via 
El X El allowed transitions. 

o Both longitudinal (Ti) and transverse (T2, T3) relaxation times must be longer than 
the PSR characteristic time scale, the time required to develop collective and coherent 
de-excitation process. 

o The initial atomic state must be prepared to have large coherence. 

o The atomic density in the initial state must exceed a certain threshold value. 

At present we plan to realize a PSR experiment with para- hydrogen (PH2) molecules, and 
details of the experiment will be described in the subsequent section. We would like to end 
this topic by emphasizing two characteristic features of the PSR experiment. First, it is a type 
of an experiment in which one tries to realize required experimental parameters or conditions. 
This is because, due to its highly non-linear nature, PSR event has a clear threshold above 
which signal can be observed without fail. This is in contrast to an experiment in which 
events occur proportional to elapse time. Second, the role of trigger laser is important as 
in the case of RENP. As stated already, the trigger laser initiates collective and coherent 
de-excitation process. Without the trigger laser, PSR events would not be initiated; this is in 
sharp contrast to the Dicke type super-radiance, in which quantum fluctuation can actually 
start coherent de-excitation process. 

4.1.3. Experimental strategy toward neutrino spectroscopy with atoms . There are two 
crucial features in realizing the neutrino spectroscopy with atoms/molecules. One is the 
experimental proof of the macro-coherent amplification mechanism, and the other is forma- 
tion of the stable condensate; both are related to the dynamical RENP factor r]^{t). Thus 
our strategy is to establish macro-coherent amplification mechanism by PSR first. As was 
stated already, we plan to realize PSR experiments with vibrational levels of pH2 molecules. 
The current status of PSR experiment with pH2 is described in detail in Sec. 4.2. 

Studies on condensate formation are one of the crucial step after establishing the PSR 
process. The points of condensate studies are 
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Fig. 26 Energy levels of pH2 relevant to the PSR experiment. 



o condensate formation method 

o time development of stored energy and coherence 

o experimental stability against PSR. 

These studies will be performed with Xe atom and other heavy atoms/molecules. Simulta- 
neously with this condensate study, target selection studies for the neutrino spectroscopy 
will be pursued. The requirements for the target are qualitatively the same as in the case of 
PSR except that the intermediate state(s) must connect the initial and final states via El 
and Ml transitions. In addition, the initial and final energy difference has profound impact 
on the neutrino spectroscopy since important physics objectives such as Dirac-Majorana 
distinction and/or CP violating phase determination can only be achieved with energy dif- 
ference comparable to the neutrino mass scale (< 0.5 eV). Experimental studies on those 
new challenging subjects are described in Sec. 4.3. Results of the some preliminary studies 
using other methods are presented in Appendix E. 



4-2. PSR experiment with para-hydrogen molecule 

In this Section we describe PSR experiment with pll2 in detail. The diagram of the pll2 
energy levels relevant to this experiment is shown in Fig. 26. The two-photon process from 
the vibrationally excited state {Xv = 1) to the ground state (Xv = 0) is measured. We note 
vibrational transitions of the El type is strictly forbidden in homo-nuclear diatomic system. 
Figure 27 shows a schematic diagram of the experimental setup. It mainly consists of three 
parts: (a) the target, (b) the excitation laser system, and (c) the trigger laser system. The 
target is LN2-temperature pll2 gas contained in a cell. The initial v = 1 state is prepared by 
the excitation laser system. The system consists of two color pulse lasers whose wavelengths 
are 532 nm and 683 nm. The trigger laser is a continuous wave (CW) laser and its wavelength 
is 4.8 nm. The excitation laser beams are combined with the trigger laser beam by a dichroic 
mirror (DCM) so that they propagate collinearly in the target. PSR outputs (left moving 
outputs) and excitation lasers are dispersed by a prism and separately detected with a photo- 
diode (PD). The backward components (i.e. right moving outputs) pass through the Faraday 
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Fig. 27 Schematic diagram of the PSR experiment with pH2. FR: Faraday rotator, A/2: 
half-wave plate, PD: photo-diode, M: mirror, DCM: dichroic mirrors, PBS: polarizing beam 
splitter. 



rotator and are reflected by the polarizing beam splitter; then they are detected by another 
PD. In the following, details of the target and laser systems are described. 

4-2.1. Gas and solid para-hydrogen target. Preparation of target, namely selection of 
atom/molecule species, choice of transition levels, preparation of initial conditions etc., is 
decisive factor for successful PSR experiments. From the point of view of atom/molecule 
selection, important parameters to be considered are: 

o Number density of atoms/molecules 
o Target length 

o Damage threshold (the maximum laser intensity which a medium can receive without 

damage) 
o Relaxation time ( Ti , T2). 

Particularly important are the number of excited atoms/molecules and coherence of |e) and 
\g). Both are determined by the number density of ground state atoms/molecules, target 
length and the excitation laser power. A higher density, a longer target and a higher dam- 
age threshold provide better conditions. However, the high-density and long- coherence are 
contradictory, because frequent interaction in high-density samples makes the coherent time 
short; thus compromise is essential for this purpose. 

For PSR, relevant transition between two levels should be El forbidden and El x El 
allowed to have a large two-photon transition moment. Vibrational transitions of homo- 
nuclear diatomic molecules fulfill these conditions. Numerical simulation in Sec. 2 and Sec. 
4.1 suggests that vibrational transition of para-hydrogen (pH2) with a large initial coherence 
is one of the ideal target for PSR experiments. Below properties of both gaseous and solid 
pH2 are described and compared. Atoms or molecules trapped in solids (pH2 for example) are 
also interesting targets for PSR/RENP; experimental studies in this direction are presented 
in Appendix E. 
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Fig. 28 Ortho-para ratio (o/p) of ther- 
mally equilibrium hydrogen at to 300 
K. The o/p can be controlled from almost 
zero to 0.75 by magnetic catalysts kept at 
corresponding temperature. 
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Fig. 29 Phase diagram of eH2.[26] Tem- 
perature and pressure at the triple point of 
eH2 is about 13.8 K and 7000 Pa, respec- 
tively. Structure of solid pH2 is h.c.p. at low- 
temperature, low o/p. 



Properties of molecular hydrogen. The molecular hydrogen, the simplest diatomic 
molecules, has attracted great interest in various fields, both theoretical and experimen- 
tal. Physical properties of the molecular hydrogen are well summarized in P. Clark Souers' 
book. [26] The nuclear spin isomers of hydrogen molecules show different nature. The nuclear 
spin of proton (/= 1/2) causes two spin states of molecular hydrogen, para (/= 0) and ortho 
(/ =1). The nuclear spin wavefunction of pH2 is anti-symmetric under the exchange of two 
identical protons while that of ortho-hydrogen (0H2) is symmetric. Therefore, the rotational 
quantum number J of pH2 and 0H2 should be even and odd in the vibronic ground state, 
respectively, because overall wavefunction should be anti-symmetric under the exchange of 
nuclei. Conversion between these two states is very slow without external magnetic fields so 
that it is possible to treat them as different molecules. On the other hand, ortho-para ratio 
(o/p) easily reaches thermal equilibrium with a magnetic catalyst. The o/p ratio of ther- 
mally equilibrium hydrogen (denoted by eH2) at room temperature is 3:1 and it converges to 
zero at the limit of K (Fig. 28) and can be controlled from almost zero to 0.75 by passing 
gaseous or liquid hydrogen into magnetic catalysts kept at corresponding temperature. It is 
worth noting that even at the triple point (13.8 K) o/p 10^^, which is the lower limit by 
this method. 

The phase diagram of eH2 is shown in Fig. 29. The temperature and pressure at the triple 
point of eH2 is about 13.8 K and 7000 Pa, respectively. The low melting point allows us to 
use gaseous sample at wide temperature ranges while solid sample should be kept around 
liquid Helium temperature. Structure of solid pH2 is known to be hexagonal closed packing 
(h.c.p.) at low temperature and low o/p. Spectroscopic data of molecular hydrogen are 
also well known. The light mass makes the rotational constants and vibrational frequency 
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fairly large. Energy splitting between J = and J = 2 (pH2) and J = 1 and J =3 (0H2) 
in the vibronic ground states is 354.38 cm~^ (510 K) and 587.06 cm~^ (845 K). This large 
energy separation strongly suppresses the population of rotationally excited states in thermal 
equilibrium at low temperatures. Even at 77 K, the Boltzman factor of the J = 2 state is 
less than 10~^. The vibrational frequency is 4401.2 cm~^. The first electronic excited state 
lies at 91700 cm~^ above the ground state and can be treated as a far-rcsonant state in our 
experimental scheme. Rovibronic energy states are not largely changed in solid pH2 because 
of weak inter-molecular interaction. Hydrogen is homo-nuclear diatomic molecule with no 
permanent electric dipole moment so that both pure rotational transition and rovibrational 
transition arc dipole-forbidden although weak electric quadrupole transition between v = 1 
and V = was observed with long optical path. [27] On the other hand, Raman transition is 
allowed and has been studied extensively, especially vibrational Raman transition. Because 
of the large vibrational energy, molecular hydrogen is often used as a medium of a Raman 
shifter. Both solid and gaseous hydrogen have been investigated as coherent Raman media by 
adiabatic manipulations. [28] [29] [30] [31] These approaches are suggestive because coherence 
of system is constitutive for PSR. 

Gaseous hydrogen. The coherence time of vibrational states of gaseous hydrogen can be 
estimated from the linewidth of vibrational Raman transitions. The linewidth depends on 
the number density (pressure) and temperature. Hereafter, the relation of number density 
and coherent time is discussed. 
In general, linewidth of gaseous sample is determined by 

o Natural lifetime 

o Doppler broadening 

o Collisions (Pressure broadening) 

o Power broadening 

o Transit-time broadening. 

The power broadening and transit-time broadening is not intrinsic linewidth of sample. It 
is to be noted that the type and nature of decoherence (homogeneous or inhomogeneous) 
are not distinguished here and only total coherence time is discussed. The parameter Ti 
describes the lifetime due to the spontaneous emission and represented by the Einstein's A 
coefficient of samples. In the case of molecular hydrogen gas, its rovibrational transition is 
dipole-forbidden. It means that the lifetime of rovibrational excited states is substantially 
long. Actually, lifetime due to the quadrupole transition of the v = 1 state is in the order of 
10® sec. [27] Natural lifetime broadening can be neglected in what follows. 

Doppler broadening, which is attributed to inhomogeneous Doppler effect caused by 
Maxwell-Boltzman distribution of molecular velocity, is dominant in low pressure gases. 
Neglecting a small relativistic correction, the full width at half maximum (FWHM) of 
Doppler broadening is given by 



where vq, T, m and c is the transition frequency, temperature, mass of molecules and the 
velocity of light, respectively. 




(80) 
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Fig. 30 Linewidth of Qi(0) Raman transition of gaseous pH2 at 81 K with ortho-para 
ratio of 1:7.7 and 3:1 function of density of pH2. 

Colhsional broadening, also called as pressure broadening, is dominant in high pressure 
gas. Detailed analysis of their lineshape is difficult because complicated inter-molecular 
interactions should be considered. [32] However, its lineshape is approximately Lorentzian 
and its FWHM is proportional to the pressure. Pressure-broadening coefficients are 
species-dependent values and typically 10 MHz per Torr. 

In the case of rovibrational transitions of molecular hydrogen gas, the pressure and tem- 
perature dependence of Raman linewidth is well studied. Bischel and Dyer [33] reported 
experimental linewidth and fitting coefficients using phenomenological equation as 



where p is the total number density of pH2 and 0H2 in amagat (1 amagat = 2.69 x 10^^ cm~^) 
and A, B is fitting coefficients. (Table 3). The second term corresponds to pressure broad- 
ening. The ffist term represent effect of Doppler broadening and " Dicke-narrowing" . The 
Dicke-narrowing appears between the Doppler dominant region and the collision dominant 
region. This narrowing is due to a smaller Doppler shift caused by frequent velocity-changing 
collisions. Fig. 30 shows linewidth of the Qi(0) Raman transition at 81 K and ortho:para = 
1:7.7 and 3:1 as a function of density of pH2. The minimum linewidth of the o/p = 1:7.7 
sample at 81 K is 117.4 MHz (corresponds to 2.7 ns) at 1.15 amagat (= 3.09 x 10^^ cm~'^). 
It is easy to prepare ~ 10^^ cm~^ pH2 gas. However, T2 is in the order of 0.1 ns at such high 
pressure. 

Damage threshold by laser field of gaseous hydrogen depends on frequency of the applied 
field and is probably determined by non-resonant multi-photon ionization in most of our 
frequency region. Recent ultrafast experiments reported multi-photon ionization of hydrogen 
molecules by visible light is dominant at > ~ TW cm~^.[34] It is easy to elongate the sample 
cell of gaseous hydrogen to 1 m or more and possible to achieve longer pass length by using 
a multi-pass cell or a cavity. 
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Table 3 Linewidth coefHcients of Raman Qi(0) transition of gaseous pH2.[33] 



T(K) 


ortho:para 


A (MHz amagat) 


B (MHz /amagat) 


81 


3:1 


189 ± 40 


29.0 ± 1.0 


81 


1:7.7 


76 ± 6 


45.4 ± 0.8 


298 


3:1 


257 lb 12 


76.6 lb 0.8 



Table 4 Typical properties of gaseous and solid pH2 Parameters of gas is at the minimum 
linewidth condition at 81 K. Parameters of the solid is at around 4 K. 

gas solid 
Density (cm-^) lO^'^ ~ 10' 

Ti (s) 10^ 
T2 (ns) 3 ~ 0.1 

Damage threshold (W cm~^) ~ 10^^ 
Target length (m) 0.1 '-^ 1 



2.6 X 10^ 
4 X 10-5 
~10 
~ 108 
< 0.1 



Solid hydrogen. Solid pH2 is an attractive target from the point of coherent experiments 
because it fulfills high density and long coherence simultaneously. The number density of 
saturated solid pH2 is about 2.6 x 10^^ cm'^ at 4 K, which corresponds to that of gaseous 
sample at 1000 atm, 300 K. Due to weak interaction, not only vibrational motion but also 
rotational motion of hydrogen are quantized and coherence time is much longer than classical 
solids. The long coherence time of the excited vibrational state is estimated to be on the order 
of 10 ns from linewidth of stimulated Raman spectroscopy. [35] The time-resolved coherent 
anti-Stokes Raman spectroscopy (TRCARS) also supported this value. Introduction of the 
TRCARS experiment and discussion of the long T2 of solid pH2 are in Appendix D. Long 
Ti (not radiative) of the v = 1 state in solid pH2 was also reported to be ~ 40 jis at 4.8 
K.[36] 

However, damage threshold of solid pH2 was reported to be 180 MW cm""^, which is 
fourth order smaller than that of gaseous sample. [30] Furthermore, it may be troublesome 
to prepare longer solid pH2 than 10 cm with optically transparent quality. Application of 
multi-pass or cavity is also difficult because of scattering in the solid. 

Comparison between pH2 gas and solid targets. Table 4 lists a set of parameters and their 
typical values relevant to the PSR experiment, comparing between gaseous and solid pH2. 
As seen, the solid phase is better from the view points of the number density and de-phasing 
time r2. Demerit of using solid phase is its low damage threshold. It limits attainable number 
density and initial coherence low; actually too low to observe PSR events with our current 
technique. On the other hand, the numerical simulation of Sec. 4.1 shows that the linear 
regime PSR may well be observed with gas phase pH2. We have thus chosen gas phase pH2 
aiming at the first observation of PSR events. It should be noted however that solid pH2 is 
much more attractive once the damage threshold limitation is overcome. Some development 
efforts along this line are described below and Appendix D. 
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Fig. 31 The cryostat for gas experiments (Infrared Laboratories, ltd. HDL-8). Both liquid 
N2 (77 K) and He (4.2 K) can be used as cryogenic liquids. Temperature of a cold stage can 
be controlled by a heater with a PID circuit. There are two view ports with IR windows for 
optical experiments. 



Experimental techniques. In solid pH2, the coherence time depends largely on the 0H2 
concentration. The linewidth of vibrational Raman transition to the v =1 state is about 
10 MHz at o/p of 2000 ppm while this becomes 60 MHz at 20000 ppm.[35] Then, a highly 
purified pH2 is desired. Purity of pH2 is also important in gaseous pH2 because FWHM of 
gaseous hydrogen at the same pH2 density is probably smaller for the pure pH2 sample as 
seen in Fig. 30, although there is no precise data at various ortho-para ratios. Almost pure 
pH2 can be obtained by passing hydrogen gas into magnetic catalysts such as FeO(OH) kept 
at about 14 K. The o/p ratio of prepared samples can be estimated in several ways. By 
measuring IR absorption spectra of solid hydrogen, the o/p ratio is determined by intensity 
of 0H2 induced pH2 vibrational transition, which can be observed in solid phase. A typical 
value of ortho-para ratio in our laboratory is about 500 ppm. 

For gaseous experiments, a sample cell was attached to a cold stage of a cryostat (Fig. 31). 
The temperature of cell was measured by Si diodes and controlled by a heater with a PID 
circuit. The cell was made of oxygen- free copper and both sides were sealed by IR windows 
with indium gaskets for an optical path. We usually use BaF2 windows for both cell and 
cryostat because of its excellent IR transmittance. The optical path is 7.5 cm in length and 
2 cm in diameter (Fig. 32). Sample gas was installed via stainless tubing connected to cell 
with silver soldering. Pressure capacity of the cell is around 500 kPa. For higher pressure 
experiments, it is better to use small diameter cell and CaF2 or Sapphire windows. For 
example, pressure capacity of a cell with 4 mm diameter and windows of 2 mm thickness is 
8 MPa (CaF2 windows) or 40 MPa (Sapphire windows). 
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Fig. 32 The home-made gas cell of 7.5 cm length and 2 cm diameter applicable to ~ 1 
K and ~ 500 kPa. The cell was made of oxygen-free copper and both sides were sealed by 
IR windows with indium gaskets for an optical path. Sample gas was installed via stainless 
tubing connected to cell with silver soldering. 

For solid experiments, a vacuum chamber with a liquid Helium bath or a closed cycle 
refrigerator (cooling power should be larger than 1.0 W at 4 K) are used. The optical grade 
crystal is prepared typically by the closed cell method. This method was originally developed 
by Oka et al..[37] The prepared almost pure pH2 gas is introduced through a thin stainless 
steel tube into a copper cell cooled by a cryostat. The cell is similar to that shown in Fig. 
32. The inlet gas flow rate is controlled to keep the cell temperature at about 8 K. The 
solid pH2 grows radially inward from the cell wall and forms a polycrystalline aggregate of 
h.c.p. crystals having their c axes normal to the wall of the cell. The sample thus grown 
is transparent. The crystals prepared by this method sometimes suffer from the cracking 
especially during the cooling cycle after the crystal growth and the inhomogeneity due to 
the polycrystalline structure. 

To prepare higher quality crystals, pressurized liquid method was developed by Hakuta et 
al.[38] Crystals in the closed cell method are obtained by direct condensation from gaseous 
hydrogen. In contrast, this method produces crystals from the pressurized liquid phase by 
which crack free single crystal can be obtained. Solid pH2 on the melting curve about 15 K 
has a similar molar volume as the solid at ~ 4 K under the saturated vapor pressure. Thus, 
a homogeneous crystal is grown at about 15 K by pressurizing the liquid-phase first, and 
then the transparent crystal is cooled down to the experimental temperature. The obtained 
pH2 crystal is transparent and uncracked and have a higher damage threshold for laser 
irradiation. 

4-2.2. Laser system for coherence preparation and trigger. In this section, we explain laser 
systems for the PSR experiment using para-H2 vibrational states. We explain the excitation 
laser system first, and the trigger laser system is described later. Since we employ a nonlinear 
optical frequency-conversion technique for the laser systems, we also explain this technique 
briefly. For detail, the reader is referred to a standard textbook such as Ref. [39] or Ref. [40]. 

As described in Sec. 4.1, it is important to prepare a large atomic (or molecular) coherence 
in the ground state | g ) and the excited state | e ) for the experiment. To prepare a large 
coherence, we manipulate the stimulated Raman process adiabatically as depicted in Fig. 
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Fig. 33 Energy level diagram to prepare coherence peg in an atomic or molecular system. 
Two color driving laser fields at frequencies of and ojc are applied. A small detuning 
of 5 is negative in this diagram. When the detuning 5 is appropriate and the intensity of 
driving lasers is sufficiently large, the magnitude of the coherence of this transition becomes 
its maximum value. 

33. For adiabatic manipulation, we apply two-color strong driving lasers whose frequencies 
are slightly detuned from the two-photon resonance of a selected Raman transition. When 
the detuning 5 is appropriate and the intensity of driving lasers is sufficiently large, the 
magnitude of the coherence of this transition becomes its maximum value \peg\ =0.5. 

In the experiment, high power and narrow-linewidth lasers are required. The time duration 
of incident laser on the order of the target relaxation time (i.e. T2) is sufficient. Thus a 
nanosecond pulse laser is preferred for two reasons. First, its power is sufficiently high: a 
peak intensity of >1 GW/cm^ with a modest beam radius of 0.1 mm can easily be obtained, 
much higher than a continuous-wave (CW) laser. Second, the Fourier transform limited 
linewidth is on the order of 100 MHz, which is sufficiently narrow for the experiment. Also 
a good temporal overlap of these two lasers is important. Figure 34 depicts a comparison of 
the beam quality and wavelength coverage range of the available technologies. The choice 
of the most relevant technology for our experiment is a crucial issue. For one of the driving 
lasers, an injection-seeded Nd:YAG laser is a good choice because of its availability of high 
power (up to several J, peak power of hundreds MW) and narrow linewidth. The wavelength 
of the other driving laser is determined by the target atom or molecule. Here, we consider 
hydrogen molecule as a target system for the PSR experiment. When we choose the first 
vibrational level as the excited state | e) (energy of 4160 cm~^) and a second harmonic of 
Nd:YAG laser (wavelength of 532 nm) as one of the driving laser, the wavelength of the 
other laser is 683 nm. 

Two types of laser techniques are commercially available at this wavelength range: a dye 
laser and optical parametric oscillator (OPO). A dye laser may be often used, however, such 
liquid laser is in general unstable and its linewidth is very broad. Usually the linewidth of 
a dye laser is on the order of GHz: 20-30 times broader than the transform limit. Another 
disadvantage of the dye laser is that frequent maintenance (replacing the dye solution and re- 
optimization of laser cavity) is required. An OPO is solid-state system and it offers very wide 
tunability (205 nm to 2550 nm). However the linewidth of such OPO is much broader than 
that of dye laser. Typically it is about 75 GHz: 300-500 times broader than the transform 
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Fig. 34 Comparison of the beam quality (linewidth and temporal profile) and the spectral 
coverage range of available laser technique. OPO: optical parametric oscillator; OPG: optical 
parametric generator; Ti:S: titanium sapphire laser; Nd:YAG: neodymium-doped yttrium 
aluminum garnet laser. The fundamental wavelength is 1064 nm and its second harmonic 
(SH) is 532 nm; FT limit: Fourier-transform limit. The wavelength tuning range of dye lasers 
and OPOs are very wide, however, its linewidth is very broad: typically 20-500 times broader 
than the FT-limit. Injection seeding technique is useful to reduce linewidths of those light 
source. 



limit. Alternatively, a titanium sapphire (Ti:S) laser may be another choice. It is a widely 
tunable all-solid-state laser system. In particular, injection seeded Ti:S laser offers Fourier 
transform limit linewidth and good beam quality. Dual-wavelength injection seeded Ti:S 
laser was used to produce near maximum molecular coherence [41]. In this case, however, 
rotational Raman transition was used. Lasing wavelength range of Ti:S laser is limited from 
700 nm to 900 nm (about 3000 cm~^). Thus this laser may be difficult to use for vibrational 
Raman process of pH2 because the Raman transition frequency is 4160 cm~^. 

Whereas those light sources mentioned above are unsuitable for the experiment, injection 
seeded OPO or optical parametric generator (OPG) is highly suitable technique. Continuous- 
wave (CW) laser seeding technique is often used to generate narrow linewidth coherent light 
[42-51]. Injection seeding requires an external light source of modest power (from a few 
I^W to a few mW) and good spectral properties; also the tuning bandwidth is limited by 
seed laser. However it offers the advantage of transform limited or nearly transform limited 
linewidth. 

An injection seeding OPG uses the process of optical parametric amplification, which is 
shown in Fig. 35 (a). A strong pump wave at the frequency of ujp and a weak signal wave at 
the frequency of cOg interact in a nonlinear optical medium. As a result, an output wave at 
the difference frequency uji = ojp — ajg is generated, and the signal wave is amplified. When 
the nonlinear medium is placed in an optical cavity, the presence of feedback produces 
oscillation. In this case, the system is called OPO. Experimentally, OPGs are especially 
preferable for injection seeding because an OPG requires no optical cavity. Thus there is no 
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Fig. 35 (a) Schematic diagram of the injection-seeded optical parametric generation. 
Corresponding energy level diagram is also shown. Schematic representations of a parametric 
frequency conversion with (b) an angle-phase-matching homogeneous bulk crystal and (c) 
a periodically poled material. In (b), the polarization of the signal and the idler waves are 
taken as extraordinary. The generated signal and idler waves at each point in the crystal does 
not propagate coUinearly with each other. In (c), the optical axis of the crystal alternates in 
orientation with period A (on the order of 5 /im to 10 //m). This case all waves collinearly 
propagate the crystal. 



need to actively control the cavity resonant frequency to match the frequency of the seeding 
laser. This makes the whole system simple. 

To obtain a high conversion efficiency, a phase matching condition should be fulfilled: 
kp = ks + ki, where kj are the wave vectors corresponding to the waves with frequencies LVj 
[j = •5; i)- The phase matching condition is interpreted as momentum conservation for the 
three photons participating in the nonlinear mixing process. Conventionally the condition 
is fulfilled in birefringent crystals that have two different refractive indices Hq and for 
the ordinary and the extraordinary waves. The ordinary wave is polarized perpendicular to 
the plane containing the propagation vector k and the optical axis. On the other hand, the 
extraordinary wave is polarized in the plane containing k and the optical axis. While the 
ordinary refractive index Uq does not depend on the propagation direction, the extraordinary 
index Uf, depends on the angle between the optical axis and k. Therefore careful tuning 
of the angle (i.e. the refractive index) is required to establish the phase matching condition. 
A serious drawback of using angle tuning is walkoff effect, which is often observed in a 
birefringent crystal[52]. Whenever the angle 9 is different from or 90 degrees, the Poynting 
vector S and the propagation vector k are not parallel for extraordinary waves. As a result, 
ordinary and extraordinary waves with parallel propagation vectors quickly diverge from 
one another as they propagate through the crystal, as depicted in Fig. 35 (b). This walkoff 
effect limits the spatial overlap of the two waves and decreases the efficiency. Thus in general 
one cannot achieve a large single-pass parametric gain with an angle-phase-matching bulk 
crystal even in the perfect phase-matching condition. To overcome low conversion efficiency, 
many systems that have been reported in literature used optical cavity to increase the pump 
power and the effective interaction length. Also a pair of crystals is often used in order to 
compensate walkoff. 
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There is a relatively new technique known as quasi-phase-matching. In quasi-pliase- 
matching, a periodically-poled material whose optical axis is inverted periodically as a 
function of position is used. The idea of quasi-phase-matching is schematically illustrated in 
Fig. 35 (c). The advantage of quasi- phase-matching is that there is no walkofF effect. Thus 
one can use a longer crystal than an angle phase-matching bulk crystal and therefore can get 
higher single-pass conversion efficiency. Since the single-pass conversion efficiency is large, 
an optical cavity is not necessarily required. Therefore an injection seeding OPG with a 
periodically poled crystal offers great simplicity of the system. 

A drawback of periodically poled crystals, as compared to conventional bulk crystals, is 
that it is difficult to make a thick crystal with present technology. This is because a high 
electric field (above 20 kV/mm) must be applied to the crystal during the fabrication process 
of periodical inversion of crystal polarization. There is no such restriction for bulk crystals, 
thus fabrication of thick crystals is feasible. Since the optical parametric process uses a 
nonlinear effect, high conversion efficiency can be obtained with high input pump intensity. 
However, there is a certain damage threshold that limits maximum input intensity for all 
nonlinear materials. In order to keep input intensity less than the damage threshold, a thicker 
crystal (i.e. larger aperture crystal) is much better for high power operation. 

There are various types of periodically poled crystals such as lithium niobate (LiNbOs, 
LN), lithium tantalate (LiTaOa, LT), and potassium titanyl phosphate (KTiOP04, KTP). 
Among these crystals, LN has the largest nonlinear coefficient. However it also has a lower 
photo-induced damage resistivity. Doping of MgO improves resistivity, but its damage 
threshold is still not very high. The nonlinear coefficient of KTP is lower than that of LN, 
and the damage threshold is almost the same as that of LN. On the other hand, stoichio- 
metric lithium tantalate (SLT) has a higher damage resistivity than LN and KTP. Although 
its nonlinear coefficient is about 2/3 of LN, the damage threshold is twice as large as LN 
[53]. Thus one can expect a higher output power by using SLT. 

We employed MgO doped PPSLT crystal (Oxide Corp., Q1532-O001) for our OPG system. 
The thickness of PPSLT is limited to 1 mm, so we used a 8 mm wide crystal in order to reduce 
the input intensity. A schematic diagram and photograph of the OPG system is shown in 
Fig. 36 (a) and (b). The PPSLT crystal was pumped by the second harmonic of a Q-switched 
injection-seeded Nd:YAG laser (Litron LPY642T). The Nd:YAG laser is flash lamp pumped 
at repetition rate of 10 Hz with a true TEMoo single transverse mode (M^ < 1.3) and single 
longitudinal mode output. Pulse energies up to 130 mJ with a pulse duration of 8 ns is 
available at 532 nm. The pump beam is expanded to elliptical shape (0.68 x 3.7 mm^) by 
means of cylindrical lenses. For injection seeding, an extended cavity diode laser (ECDL) 
in the Littrow configuration was used. We used a commercially available laser diode chip 
(TOPTICA LD-0685-0050-3, no anti-reflection (AR) coating) for the ECDL. The output 
power of the ECDL is more than 10 mW with a typical mode-hop-free scanning range 
of 3 GHz. A polarization maintaining optical fiber was used for transverse mode cleaning. 
Typically, a fiber output of 2 mW was used for injection seeding. The injection seeding beam 
was also expanded to elliptical shape (1x8 mm^). The pump and injection seeding lasers 
are combined with a dichroic mirror. The PPSLT crystal is 40 mm long x 1 mm thick x 
8 mm wide (periodic domain width of 7 mm) and periodicity of 10.3 /xm. Phase matching 
temperature at 683 nm was measured to be around 110°C. The crystal is AR coated for 
three wavelengths: pump (532 nm), signal (683 nm), and idler (2.4 /xm). Figure 37 shows a 
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Fig. 36 (a) Schematic diagram and (b) Photograph of the excitation laser system. ECDL: 
extended cavity laser diode; DL: diode laser chip; Nd:YAG: pulsed Nd:YAG laser; L: spherical 
lens; CLs: cylindrical lenses; Dl: dichroic mirror used to combine pump laser (Nd:YAG) and 
continuous- wave seed laser (ECDL); D2-D3: dichroic mirrors used to separate the pump plus 
the signal beam and the idler (2.4 /im); BD: beam dumper; Iso: isolator; PMF: polarization 
maintaining optical fiber; PPSLT: periodically-poled stoichiometric lithium tantalate crystal. 
The PPSLT crystal is pumped by the second harmonic of a Q-switched injection-seeded 
Nd:YAG laser, and generates signal (683 nm) and idler waves. 



typical temporal profile of the transmitted pump beam and generated signal (683 nm) light. 
Good temporal overlap is obtained. 

We next describe the trigger laser system. According to the PSR simulation (Sec. 4.1), 
the power of trigger laser is less important, however the narrow linewidth (<MHz) may 
be crucial. The wavelength of the trigger laser is 4.8 fiui. While there are many lasers 
available at visible or near-infrared (NIR) region, coherent light sources emitting mid-infrared 
(MIR) radiation is less common. Recently, quantum cascade lasers, which directly radiate 
mid-infrared wavelength light, are commercially available. However the linewidth is broad 
(about 45 MHz) and they are still expensive. Moreover, peripheral optical devices in the 
MIR spectral region such as optical isolators, gratings, high-reflectivity low-loss mirrors for 
high-finesse optical cavity, or a wavelength meter are also less common. Thus there are 
practical drawbacks when one develops laser sources with direct laser radiation devices in 
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Fig. 37 Temporal profile of OPG output. Green line shows the temporal profile of the 
pump laser transmitted the crystal. Red line shows that of the generated signal (683 nm) 
light. 

the MIR region. Instead, a light source based on nonlinear optical process (i.e. difference- 
frequency generation) with a quasi-phase-matching material pumped by NIR lasers is highly 
competitive. It allows to transfer the high performance characteristics of the pump laser at 
NIR region to the MIR: precise wavelength resolution and narrow linewidth. 

Figure 38 shows the block diagram and the photograph of the trigger laser system. We 
employ two of master oscillator power amplifier (MOPA) systems as pump sources of DFG. 
The wavelength of the MIR light is determined by measuring wavelengths of both pump lasers 
by a wavelength meter (HighFinesse, WS6-200). One of the pump laser is combined with the 
other by DCM and coupled into optical fiber. We employ waveguide PP-MgO:LN crystal 
(NTT electronics, WD4800-000-A-B-C). A great advantage of waveguide PPLN device is 
that one can obtain longer interaction length with higher beam intensity, which yields two 
or three orders of magnitude larger single-pass frequency-conversion efficiencies than no- 
waveguide crystals. This is especially effective for CW laser because its beam intensity is 
relatively low. 

4-3. Towards RENP experiment with Xe 

4-3.1. RENP with Xe atoms. As described in Sec. 1, Xe atom is one of the important 
candidates for studying the RENP. The lowest excited state 5^^(^^3/2)65 ^[3/2]2 with total 
electron angular momentum J = 2 in Xe atom can decay only through M2 transition to 
the ground state 5p^(^5o), and thus this J = 2 state is metastable with lifetime of 0(40) 
s [54] . The RENP process in Xe atom would be recognized by observing a characteristic decay 
spectrum of emitted photon from the metastable state, which has distinctive threshold at the 
photon energy of 4.16 eV, half of the energy difference between the metastable and ground 
states. Figure 39 shows such typical RENP spectrum, where different threshold locations ujij 
appear depending on mass scale of participating neutrino pair. 

In order to successfully observe the RENP process, understanding and precise control of a 
coherent state of atoms interacting with static field condensate is inevitably important [8]. 
Therefore the investigation on feasibility and further research/development is needed for 
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Fig. 38 Trigger laser system: (a) Block diagram of the trigger laser system. The wave- 
length of the trigger laser is 4.8 ^m. We employ difference-frequency generation (DFG) to 
obtain coherent light at this mid-infrared wavelength range. For pump lasers of DFG, we 
employ two master oscillator power amplifier (MOPA) systems at wavelength of 871 nm and 
1064 nm, respectively. Both lasers are coupled into optical fiber and then coupled to the 
waveguide PPLN crystal. Typical output power from each device is indicated in the figure, 
(b) Photograph of the system. Note that MOPA systems are only shown in the photograph. 
ECDL: extended cavity laser diode; TA: tapered amplifier; Iso: isolator; DCM: dichroic mir- 
ror used to combine two pump lasers; FC: fiber coupler; WG-PPLN: waveguide PP-MgO:LN 
crystal; L: collimator lens; Ge: Germanium filter. 



the preparation and confirmation of the atomic coherence between the metastable and the 
ground state in Xe atom. The fundamental characteristics of metastable state in Xe, numer- 
ical calculation of achievable coherence, and the prospect for experiments are described 
below. 

4-3.2. Properties of Xe atom in gas phase. The atomic coherence between the metastable 
and the ground state in Xe can be produced by applying two light-fields ("pumping" 
and "coupling or sometime called stokes") on so-called A- type three- level atomic system. 
However, the achievable atomic coherence tends to be limited from several origins, which 
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Fig. 39 Examples of Xe RENP spectra near threshold. Assumed are Dirac neutrinos with 
the minimum masses tuq of 2 (black) , 20 (red) and 50 (blue) in the normal (solid) or inverted 
(dashed) hierarchy. 



broaden in effect the line widths in the relevant atomic transitions. Because these line widths 
are affected by physical and spectroscopic properties of Xe atom, these parameters are 
summarized here. 

The metastable state of Xe atom has long lifetime over 40 s for even isotopes in low- 
pressure gas phase. The radiative lifetimes of 42.9(9) s and 42.4(13) s for ^^^Xe and ^^^Xe 
respectively, are reported where the observation was performed in MOT (Magneto-optical 
trap) experiment [54]. The density of Xe atom can be increased up to 2.5 x 10^^ /cm^ with 
a pressure of 1 atm at room temperature. If the gas is cooled to 190 K (boiling point at 
pressure of 3 atom) and the pressure is increased to 3 atom, the number density reaches 
1.1 X lO^'^ /cm^. The atomic collision should be considered for relaxation of metastable 
atoms in relatively high density gas phase we here consider for our purpose (> 10^9 /cm^). 
The collisional relaxations (atom-atom and atom-wall) of the metastable Xe in a gas cell are 
also reported, where the observed relaxation rate takes minimum (~ 10^ s~^; depending on 
cell size) at a pressure of 0.2 - 1 Torr and reaches to 1 x 10^ s~^ at 10 Torr [55]. From this 
reported result, the collisional relaxation time of the metastable state of 0(0.1 ~ 0.01) ^s is 
expected at a pressure of 1 atm that is enough long for experimental RENP studies. 

Populating Xe atoms into the metastable state is performed through two-photon excitation 
from the ground state to 5p^{'^P^/2)^P states as shown in Fig. 40. Among these 6p-blanches 
the two states, 5/(2^3/3)6^ ^[3/2]2 and 5/(2^3/3)6^ ^[5/2]2, are allowed both for the two- 
photon excitation from the ground state and single-photon El transition to the metastable 
state in the A- type three-state system. Both of these 6p-states have El-transitions to the 
6s [3/2] 2 (metastable) and 6s [3/2] 1 which promptly decays into the ground state with El 
transition. The consecutive El decays through 6s [3/2] 1 state is one of sources for relax- 
ations and decoherences in the three-level system. The spectroscopic parameters for these 
two states including the wavelength Xtp for two-photon and Agsj-ep for the 6p-6s transi- 
tions, transition probabilities A^gj^-Qp, are summarized in Table 5. Among these two states. 
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Fig. 40 Schematic energy level for populating Xe atom into metastable state. 



Table 5 The parameters of relevant excited 6p states for populating Xe atoms into the 
metastable state 5p^(2P3/2)6s ^[3/2] 2. The values of transition probabilities ^gsj^ep are 
taken from [56]. 





(eV) 


(nm) 


^6s2-6p 
(nm) 


^6s2-(-6p 

(s-1) 


A6sl-6p 

(nm) 


(s-1) 


5/(2P3/2)6p 2[3/2]2 
5p5(2P3/2)6p 2[5/2]2 


9.821 
9.686 


252.5 
256.0 


823.4 
904.8 


2.24 X 10^ 
9.72 X IQS 


895.5 
992.6 


9.42 X 10*^ 
1.74 X lO'^ 



5^^(2^3/2)6^ ^[3/2]2 is preferable in terms of usability of wavelength in IR region. Therefore 
this state is used in the following estimation for atomic coherences. 

Another mechanism of the relaxation is the ionization of Xe atoms in excited 6p or 6s 
states due to absorption of third photon from the pump field. The first ionization energy of 
12.1 eV for Xe is only 2.3 eV higher than the 6p state. The ionization rates for these states 
are proportional to the light field strength. There are some experimental and theoretical 
studies for ionization of Xe atoms in 6p states. For the state 5p^ 6p[3/2]2, the observed 
ionization rate of Ti = 3.0 • I252 [s~^] is reported, where I252 is field strength of 252.5 nm in 
unit of W/cm^ [57]. Although the metastable Xe atoms can also absorb the third photon 
to be ionized, the ionization rate of Fi = 0.19 • I252 [s~^] roughly estimated by the reported 
photoionization cross section [58] is smaller than that of 6p state. 

It must be noted also on the Doppler broadenings in each transition which affect the phase 
coherence in the relevant two states. At room temperature, the Doppler widths (FWHM) of 
the 5p-6p and 6p-6s transitions are 2.6 GHz and 0.40 GHz, respectively. 

4-3.3. Coherence preparation and its measurement. Atomic coherence between the 
metastable and ground state in Xe atom is produced by two laser fields, one for two-photon 
excitation from ground |c) to the 6p state \a) (pump field with frequency ujp), another for 
coupling between 6p and the metastable state \b) (coupling field with frequency Wc)- As 
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described previously, the preparation process of coherence includes the two-photon excita- 
tion and some relaxations and is then complex. Therefore, it is important to quantitatively 
estimate how large atomic coherence can be achieved. We then perform numerical simulations 
for the coherence of metastablc Xe. 

The time evolution of populations and coherences for A-typc three-level system in the pres- 
ence of the two laser fields are estimated by the optical-Bloch equations. The optical Bloch 
equation is the homogeneous (spatially constant) version of the Maxwell-Bloch equation in 
which the light-field is taken as given. The equations for the populations paaj Pbb and Pcc in 
each state are described as 



Pa 



(Pab - Pba) + ^(Pac " Pea) " (r^d-aU) + vf^)p^ (82) 



,-0 ^ , p(rad-meta) 
Pbb = ^(Pba-Pab)-(ri'°'^^+rC°Vbb + ^^^ Paa (83) 



2 
2 

Pec = ^^(Pca-Pae) + "Y^Paa, (84) 

and similarly for the coherences, the following equations describe the time evolutions 

Pab = -(Tab - «^ab)Pab + ^^(Paa - Pbb) - ^^Peb (85) 
Pae = -(7ac - «^ac)Pae + ^^(Paa - Pec) - ^^^Pbe (86) 
Pcb = -(7eb-^(<^ab-<^ae))Pcb+^^Pea- ^Y^Pab- (87) 

Here the "slow" variables for coherences are introduced as 

Pab = Pabe-*("- Pae = Pa^e-'^^--^-^)*, = Pebe"*"-*, (88) 

in order to move into appropriate rotating frames. Jlac and Jlab = Pab-E'c/^ represent the 
Rabi frequencies for the pump and coupling field, respectively. Jlac is defined as an eff'ective 
two-photon Rabi frequency for the field at frequency of Wp; ilac = Oac-^'p/^) where Oae = 
X^i /^iiMj2/2/i(wi — Wp) and the [Xij indicates the corresponding dipole matrix elements. E'p 
and E'c are the amplitudes of the applied pump and coupling fields respectively, and these 
field amplitudes are described with corresponding intensities / of the incident lights as £^ = 
\/2ZoI = 27.5A//(W/cm2) V/cm. In these equations we include the spontaneous decay term 
p(rad all) T^^^^ mcta) gQj.j.ggpQj^(jjj^g decay ratcs from \a) to two 6s states and only to 
the metastable state \b). These decay rates take the values ri''^'^ ^^^^ = 2.63 x 10^ [59] and 
p(rad meta) _ 2 24 ^ 10'' s~^ [56]. The state |a) tends to decay into the ground state |c) with 
two El decays through 5p6s ^[3/2]i whose decay rate is estimated as F^™^^ = {1/Ta-^esi + 
l/r6si->c)~^ = 9-1 X 10^ s~^ [56, 60]. As described in the previous subsection, the states 
I a) and \b) also have decay ratcs Ta°^^ and r|j'°°^ through ionization by the pump field as 
relaxation processes. The coUisional loss of population in |6) is assumed r^°^^ = 5 x 10^ s^^ 
(we here assume high density; > 1 atm). In this calculation, the widths of Dopplcr broadening 
in each transition are included in the relaxation terms 'jij in coherence. These values are 
7ac = 1-6 X 10-'^° s~^, 7ab = 2.5 X 10^ s"-*^, and Tbc = 1-4 x 10^ s~^. The pressure broadening 
assumed to be 1 x 10^ at a pressure of 1 atm (room temperature) is also considered in 
these decoherence terms. We assume that the detunings Sgh and Sg^ in each field are larger 
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than the Hnewidths in the corresponding transitions, so that both of these transitions by the 
appHed fields are adiabatic. 

The dependence of population in the 6p state and the coherence pac between 6p and 5p 
states on the intensity of the incident pump light is firstly estimated by the above equations. 
The incident pump light is assumed to have a Gaussian envelope with a width of 5 ns and 
to have a peak intensity at i = 10 ns. The detuning (5ac of pump frequency is assumed to be 
1.2 X 10^° s"^ (2.0 GHz) so as to be larger than the Doppler broadening 1.3 GHz (HWHM). 
Figure 41 shows time evolutions of pac and paa when the pump field is applied. The effective 
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Fig. 41 Numerically calculated atomic coherence pac between ground state and 6p state 
(left panel) and the population paa (right panel) when only the two-photon field (252 nm) 
are applied. 



two- photon Rabi frequency J^ac is calculated with the effective dipole moment ai2- Here 
only 6s^[3/2]i state is taken as the intermediate state \i) for calculating the ai2- Thus riac 
is estimated as 

Oac = 5.44 • /252 s-\ (89) 

where I252 denotes the pump- field strength in unit of W/cm^. These results shows that the 
effect of ionization loss starts to appear with the field strength of I252 ~ 1 GW/cm^, and that 
the coherence achieves /9ac ~ 0.06 at that field. The achievable coherence tends to saturate 
around pac ~ 0.1 due to the ionization effect. 

The coherence pbc between the metastable and the ground state is calculated by adding 
the coupling field of 823 nm to the above situation. The Rabi frequency for this transition 
can be calculated with the ^-coefficient A = 2.24 x 10^ [56] in the transition as 

Oab = 5.43 X 10« • /fe^ s'^ (90) 

where I823 is the coupling- field strength in unit of W/cm^. The estimated atomic coherences 
Pbc with /823 = 10 ~ 5000 W/cm^ and I252 = 500 MW/cm^ are summarized in Fig. 42. The 
populations /?bb in the metastable state are also shown in this Figure. The detuning (5bc 
in coupling field is assumed to be 2.5 x 10^ s~^ (0.4 GHz) which is larger than Doppler 
broadening 0.2 GHz (HWHM). The expected coherence increases to pbc = 0.01 ~ 0.1 when 
the coupling Rabi frequency ilab exceeds the two-photon Rabi frequency ilac- As shown in 
the results, the coherence phc = 0.05 is available if the pump and coupling fields are increased 
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Fig. 42 Numerically calculated atomic coherence p\,c between the the metastable and 
ground states. 



to /250 = 500 MW/cm^ and /823 = 5 kW/cm^. In this condition each Rabi frequency takes 
^ab = 3.8 X 10^*^ and J^ac = 2.7 x 10^ s^^. The population /Obb in the metastable state 
reaches ~ 0.05 under this condition. 

There would not be straightforward way to confirm presence of coherence /?bc! because 
there is no adequate transition path between \h) and |c) for coherence measurements such 
as coherent anti-stokes Raman scattering (CARS). This is a same situation for both of the 
coherence between 6s2 and (|c)), and between 6si and 5p^. However, it is known 
that the three level system behaves as the dressed state system when the strong coupling field 
is applied (Rabi frequency exceeds the line broadening), and this strong coupling between 
these states can be confirmed by observing that absorption of the radiation from the state |6) 
to |c) is reduced, called as Electromagnetically Induced Transparency (EIT) [61]. A relevant 
case for this phenomenon was studied with atomic hydrogen [62] which has a similar A- 
type three-level structure as shown in Fig. 43. In both of these systems, the atoms are 
pumped to the excited state by two-photon transition and the atoms tend to be ionized 
by the applied fields. In the EIT experiments with hydrogen atom investigated in [62], the 
transparency of the VUV radiation (103 nm) from the 3p to the ground state and the 
reduction of the photoion signal from the atomic system were observed when the coupling 
Rabi frequency 02s-3p exceeds the line broadening ^Doppier = 2.3 x 10^^ s^^. In Xe atom, 
this method can be applied to the three- level system of |a), |c) and 6s ^[3/2]i (instead of 
1 6)), where the transparency of the VUV radiation (147 nm) from the 6s ^[3/2]i is expected 
to be observed. The required condition for the applied coupling field in this experiment is 
^^coupiing > ^Doppier ~ 2 X lO^'^ s~^ corresponding the coupling field intensity of > 2 kW/cm^. 
The coupling between the three states |a), |5) and |c), which is the target system in our 
studies, can be estimated from the observation in the above system and the numerical 
simulation. The reduction of the photoion signal with the EIT condition would also be 
useful for evaluation of the atomic coherence as the experiments with hydrogen. 

4.3.4- Prospect for experiment. As described above, the atomic coherence of the 
metastable state in Xe atom is expected to reach ~ 0.05 which enables to study the 
RENP process by the numerical simulations where the pump and the coupling fields lead 
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Fig. 43 Comparison of the A-type three-level systems between Xe and Hydrogen atom. 
The two-photon excitation scheme and the ionization are common in these species. 



adiabatic transitions. In order to perform the experiments for the coherence preparation and 
measurement, the following experimental conditions are required: 

o Linewidth of the applied fields is narrower than the line broadening 
o Smooth temporal profile for the applied fields 
o Smooth spacial profile for the applied fields. 

The main broadening in Xe atom in gas phase arises from the Doppler effect, 
and estimated 1.6 x lO^*' s~^ (2.6GHz(FWHM)) for two-photon transition and 2.5 x 
10*^ s"^ (0.4GHz(FWHM)) for couphng transition at room temperature. The pressure broad- 
ening would not be negligible in the high-pressure atomic gas, which becomes ~ 10^ s~^ at 
the pressure of 1 atm. The required linewidth is then < 100 MHz. Therefore the laser system 
with enough narrow linewidth such as an injection seeded Ti:Sapphire is required (Sec. 4. 2. 2). 
As listed above the applied fields should have smooth spacial profiles for both of the pump 
and coupling, and these fields should be overlapped in space. The available field strength 
should be estimated with the realistic beam diameter for this overlapping. If the beam diam- 
eters are assumed that c/'pump 

= 0.2 mm, ^coupling = 0.5 mm and the pulse energy is assumed 
10 mJ/pulse with a pulse width of 10 ns for both fields, the field strength can then reach 
-^pump = 600 MW/cm^ and /coupling = 100 MW/cm^ with which the coherence /Obc becomes 
> 0.05. 

The conceptual setup for the experiment is shown in Fig. 44. The high-quality pulsed 
beams for pump and coupling can be produced with Ti:S lasers injection seeded by Nd:YAG 
laser. The pumping field of 252 nm is obtained through third harmonic generation (THG) 
by using BBO crystals from the output beam of a Ti:S Laser. The detectors includes a 
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monochromator and a photomultiplier (PMT) for the light, and an electron multiplier (EMT) 
for the ionized Xe. The fundamental experiment as R&D for spectroscopy of the metastable 
Xe is now being performed with relatively broad linewidth lasers such as a pulsed Dye. It 
contains spectroscopy of transitions from the two-photon excited 6p states, investigation of 
populations in the metastable state, and lifetime of the metastable states. 



THG system 




Fig. 44 The conceptual setup for the coherence preparation and measurement of Xe. 



5. Summary and prospects 

In the present work we first explained in detail theoretical principles of the neutrino mass 
spectroscopy with atoms and molecules, which aims at determination of undetermined 
important neutrino parameters. The key idea for measurements is to amplify otherwise 
small rates of radiative emission of neutrino pair (RENP), |e) — >• \g) + J + i^iVj (with Vi, Vj 
mass eigenstates), by developing the macro-coherent medium polarization among target 
atoms/molecules strongly coupled to fields inside the medium, in order to stimulate the pro- 
cess cooperatively. The amplification is realized by trigger irradiation of two colors into a 
well prepared target state of good phase coherence between atoms in two relevant levels, |e) 
and \g). The amplification also works for two-photon emission called paired super-radiance 
(PSR), |e) — >• |(7) -|- 7 -|- 7. Detailed account of the master equation for PSR and coherence 
development for RENP is presented and results of numerical simulations on the pH2 Xv = 1 
vibrational transition have been presented both for explosive PSR events and events in the 
weaker linear regime. Our master equation includes effects of phase decoherence of medium 
polarization and decay of population difference. 

By selecting a metastable state |e) forbidden to decay to lower levels via El transition, 
one can obtain a large ElxEl two-photon PSR rate for |e) — >• \g). A good example of this 
feature is pH2 vibrational transition Xv = 1 — ?■ 0. The PSR event may occur explosively, as 
shown in Figure 4 of Section 1, if the target relaxation time T2 is larger than some number, 
for the target number density n = lO^^cm"^ this number being of order a few to several 
nano seconds. Even for smaller T2's, the enhanced output over the input trigger is expected 
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as discussed in Section 4. Discovery of these PSR events serves to prove our principle of 
macro-coherence, a necessary prerequisite for RENP experiments. 

The RENP spectral rate is given by the factorized formula r^2i/(w) = ^oI{'^)Voj{i)- The 
factor Po determines the overall rate in the unit of 1/time, and for Xe it is of order, 
1 Hz (n/1022cm-3)3 (y/lQ^ cm ). The spectral shape as a function of photon energy ui is 
given by I(oj), which serves to determine the neutrino parameters. Detailed photon spec- 
trum calculations of Xe and I2 transitions, having available energy 8.3 eV and 0.88 eV, 
respectively, have been given to explain how the parameter determination is made possi- 
ble. The dynamical factor r]oj{t) is the space integrated quantity over the entire target of 
\E~{Ri — ii?2)P/(4eegn^), the dimensionless squared product of the coherent polarization 
— iR2)/2 and the field E^~^ inside the target. 

The dynamical factor rj^iit) is highly time dependent in the phase of large PSR related 
activities. It asymptotically approaches a constant value of the stable state made of field 
condensate supported by coherent medium polarization. For RENP it is important to utilize 
this factor as large as possible. How the condensate state is described is explained and a 
sample simulation has been presented, leaving more detailed investigation to future works. 
If the amplification works as expected, the absolute neutrino mass determination along 
with distinction of mass hierarchy pattern (normal or inverted) becomes feasible. Other 
important neutrino parameters such as Majorana/Dirac distinction along with CP violating 
phases may also be determined if the initial metastable atomic state lies close to the ground 
level, preferably within a fraction of eV. 

Heavy atoms/molecules which have a large breaking of the LS coupling scheme are favored 
for RENP target candidates, since they may give a large enough amplitude of Mix El type 
for RENP between |e) and \g) . 

In the later half of the paper experimental status of our project has been reported, along 
with an overview of the entire experimental strategy. 

Our basic strategy towards precision neutrino mass spectroscopy is first to prove the macro- 
coherent amplification principle by the PSR process. To this end, we have chosen gas-phase 
para- hydrogen (PH2) get, and focus on the El forbidden transition from the vibra- 

tional excited state Xv=l to the ground state Xv=0. As is discussed in Sec. 2 and 4, PSR 
is highly non-linear phenomenon, and can be observed only when a certain set of initial 
conditions are realized. In the case of pH2, they are the number density {n), the initial 
coherence (ri or r2), and dephasing time {T2 and r3). We expect all necessary condition 
can be fulfilled in the current technology available to us at least for observing PSR events 
in the linear regime. As described in Sec. 4, all components are prepared, and ready for the 
experiments. 

It is crucial to achieve a large value of the RENP dynamical factor r/^j (t) to realize high pre- 
cision mass spectroscopy. We plan to study evolution of macroscopic coherence and formation 
of field condensates, both of which directly affects r]^{t), using Xe atom as an experimen- 
tal platform. Details of the current status on the Xe experiment can be found in Sec. 4. 3. 
Condensed-phase targets are highly desirable not only for RENP experiments but also for 
observation of explosive PSR events. We are seeking and trying to develop high density solid 
targets with long dephasing time. Some preliminary results of our efforts along this line are 
presented in Appendix D and E. 
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It is evident that we need much more detailed simulations and experimental studies on 
how best we should design this challenging experiment. Some entirely new idea is welcome 
and may further enlighten prospects towards this important objective of neutrino physics. 



Note added in proof. After the original submission of the manuscript some of us succeeded 
in deriving breathing solitons which are useful as an ideal target state for radiative neutrino 
pair emission (RENP). This shall be discussed elsewhere in a separate work [92]. 
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A. Electroweak interaction under nuclear Coulomb potential 

A.l. Electrons in atom and their electromagnetic interaction 

Electrons in an atom are described by solutions of the Dirac equation under the nuclear 
Coulomb potential [63] augmented by interaction with other electrons in the atom. Bound 
and nearly free (modified by the nuclear Coulomb field) electron states thus obtained form a 
complete set of states that may be used as an expansion basis of electron field denoted by a 
four component operator ip(,{x,t). The interaction of atomic electrons with the electromag- 
netic field A{x, t) is best described in physical terms by using the radiation gauge, as given 
for instance in standard textbooks of QED such as [63]. The radiation gauge field satisfies 
V • A = and possesses two transverse degrees of freedom representing two transverse waves 
of photon. In this gauge the static Coulomb interaction is separated from the transverse field 
interaction, hence is convenient for our purpose. 

In the present work we deal with interaction of atomic electrons with the transverse elec- 
tromagnetic field A and three massive neutrino field Vi. This way we can derive an effective 
hamiltonian of atomic electrons in the electroweak theory [11] we need for our purpose. 

In most cases of atomic electrons relativistic effects are minor and atomic electrons have 
bound energies miich less than the electron rest mass. Or equivalently, the average electron 
velocity within an atom is much less than the light velocity. This makes non-relativistic 
treatment sufficient, although in some cases relativistic corrections such as LS coupling are 
important. Under this circumstance the electron field operator can be decomposed into two 
parts; a large and a small two-component fields. The major large two-component field is 
denoted by a new operator (j)e{x,t). 

As is well known [63], the electromagnetic interaction of atomic electrons is governed by 
interaction hamiltonian (or rather the hamiltonian density to be spatially integrated), 

e4>l{x, ^)(-^V) • A{x, t)^^, t) ■ (Al) 

When this interaction hamiltonian is applied to atomic process, it gives a vertex amplitude 
for the atomic transition \i) — >■ j/), 

e(p\ {x, t) (-i V) • A{x, t)(pi {x, t) , (A2) 

where cpa is a wave function solving the nuclear Coulomb potential problem [17]. 

A further simplification arises by taking the long wavelength limit of emitted photons (the 
wavelength of photon much larger than the atomic size) , and this leads to what is called the 
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El transition in the literature [63], [17], where one can replace the vertex amplitude above 
by 

(p^j-{x,t)d ■ E{x,t)ipi{x,t) , d = ex. (A3) 

In the many-electron system in atom the dipole operator d should be replaced by a summed 
quantity over all contributing electrons at position Xn', d = ej^n^n- The wave functions ip 
should also be replaced by those of the multi-electron states [17]. 

A. 2. Quantized massive neutrino fields 

We have two possibilities to describe massive neutrinos: one is the Dirac type and another 
is the Majorana one. The neutrino and the anti- neutrino are distinct particles in the Dirac 
case, while they are identical in the Majorana case. RENP final states discriminate these 
cases due to the identical particle effect as explained in Sees. 1 and 3. The neutrino expansion 
base may be taken, effectively and conveniently, as plane wave functions, since the neutrino 
interaction with atomic nuclei is very weak. 
For the Dirac type, a quantized free neutrino field is represented as 

i^'^ix) = [Kp, X)uip, A)e-'f- + Sip, \)v{p, X)e'P-] , (A4) 

p,X 

where b{p, A) is the annihilation operator of the neutrino with a momentum p and a helicity A, 
and d^ {p, A) denotes the creation operator of the anti-neutrino. The four-component spinors 
u{p, A)e~*^^ and v{p, A)e*^^ are solutions of the Dirac equation. We employ the following 
convention for the neutrino momentum sum, 

where Ep = ^/m^ is the energy of neutrino of mass m. The Dirac field tp^ either anni- 
hilates a neutrino or creates an anti-neutrino. A free Majorana neutrino field is quantized 
as 

^^(x) = ^b{p, \)u{p, A)e-'f^ + b^{p, \)v{p, A)e^P^] , (A6) 

p,A 

where b (b^) represents the annihilation (creation) operator of the Majorana neutrino and 
satisfies tp^ = C{ijj^)'^ with C being the charge conjugation matrix. ^ We note that 
there is no distinction between neutrinos and anti-neutrinos in the Majorana case, and both 
the creation and annihilation operators of the Majorana neutrino appear in ijj^ . 

A. 3. Effective interaction of neutrinos with electrons 

The neutrino fields in the flavor eigen-basis are related to those in the mass eigen-basis by 
a unitary transformation as 

Va{x) = ^ Uai nix) , (A7) 

i 

where a = e, ij,,t specifies a neutrino flavor, i = 1, 2, 3 denotes a mass eigenstate, and Viix) = 
ip^{x) or %l)^{x) with mass mj. The lepton mixing matrix is written as a product of two 



^ This effectively reduces to two independent solutions in accordance with two-component 
formalism. The field quantization in the two-component formalism is given in [11] and [9]. 



64/85 





Charged Current 



-ij J 
Neutral Current 



Fig. Al Feynman diagram for W and Z exchange weak interaction of electron with 
neutrinos. 



unitary matrices [1]: 



where 



U = VP, 



V 



C12C13 



S12C13 



-1(5 



-S12C23 - Cl2S23Sl3e*^ C12C23 - Si2S23Sl3e'^ S23C13 



AS 



JS 



(A8) 



(A9) 



SI2S23 - Cl2C23Sl3e C12C23 - Sl2S23Sl3e C23C13 

with Cij = cos 9ij and Sij = s'm Oij. The diagonal unitary matrix P may be expressed by 

P = diag.(l,e^",e^^), (AlO) 

for Majorana neutrinos, and we rotate away the phases a and /3 for Dirac neutrinos. 

The effective interaction of neutrinos with the electron field arises from W and Z exchange 
CC and NC interaction [11], as depicted in Fig.Al, and is given by 



Vij - aij'j^je, 



(All) 



where 



U:^UeJ - Q - 2 sin^ Ow^ Sij, a,, = U:iUej - \kr (A12) 



The interaction above is written using the four component spinor notation. It has been 
shown that for non-relativistic atomic electrons the leading neutrino interaction stems from 
the axial-vector part and in the two-component notation it can be written using the spin 
operator of the form. 



(A13) 



This results in the magnetic type of transition operator on the electron side. We may thus 
use, for the spin matrix element of atomic transitions from a bound state |a) to their 
wave functions Lpa{f)^ '^hif)^ and compute the overlapping integral of the kind, 



(A14) 



In practice, it is better to use experimental data for this if they are available. 
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B. Mathematical structure of Maxwell-Bloch equation 
B.l. Maxwell-Bloch equation for SR 

We consider equations, Eq. 19 ~ 21 in the text, for the single photon SR. First rescale time, 
length and field by 

r=^, ^ = ^, e = -idTnE, (Bl) 

where Tr should have the dimension of time/length. It is further convenient to introduce 
the retarded time r_ = r — ^ and use (r_, 4) as independent variables. In terms of these the 
Maxwell-Bloch equation is converted to 

d-R = Ze, (B2) 
d-Z = -^{e*R + eR*), (B3) 
d^e = R, (B4) 

with Tr = 1/ V2TrcJ(Pn determined from consistency. 

Using the conservation law that follows, + Z^) = , we introduce an angle 

function 9{t-,^): 

R = e'^B{^) sin 0(t_, ^ , Z = B(0 cos e(r_, ^ ■ (B5) 

In order to avoid unnecessary complication we assume real field e and set the phase (p = 0. 
The Maxwell-Bloch equation then becomes 

^e = Bsme, (B6) 

e = d-e. (B7) 

The first equation is known as the sine-Gordon equation in 1-1-1 space-time dimensions. 

Scaling invariance under r_ — >■ /3t- , ^ , exists in this system. The Burnham-Chiao 

equation [64] is obtained using the variable, w = 2^Jt-^ = 1^ (r — ^)^, 

dH \ dd ^ . ^ 

—^ + -— = Bsme. B8 

aw'^ w aw 

Since this is an ordinary differential equation, it is much easier to solve than 1-|-1 dimensional 
partial differential equation. This solution is illustrated in Fig.Bl, which shows interesting 
ringing structure of emitted pulses. 

B.2. Maxwell-Bloch equation for PSR 

We recapitulate the master equation for the most general case of two-color trigger irradiation 

of frequencies ui + uj2 = ^eg, as given in [14]. 
We introduce the dimensionless unit: 

|^{l),{2)|2 ^ 

{^,t) = ifJ-x , iJ.t) , n{uj) = Ceguageiuj, e^g - w) , je^^^^^^^^ 1^ = , ri = — , (B9) 

„, , , ^^ _ dpedpg{epg + Cpe) 

age[iO, €eg — id) — 7 r . iBiUj 

Assume R-mover field of frequency coi and L-mover of frequency UJ2 (neither R-mover 
of frequency U2 nor L-mover of frequency uji). Note the frequency dependence of cou- 
pling parameters aab{iOi,uJ2) and the universal relation aeg{oJi,co2) = (Xge{coi,L02) for any 
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Fig. Bl Auto-modeling solution numerically computed from eq.(B8). Initial values are 
e{T = 2x 10-^) = ,9{t = 2x 10"^) = and target length = 1. Constant .B ~ 1 is 
assumed. 



combination of uJi + UJ2 = Ceg- The master equations for medium polarization and fields are 
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with Uaa defined by 



- ^2) 



(B21) 



The single mode equations in the text are readily derived by taking oji = €eg/2, Oj = 1, aij = 
1 and all 7jf a, 6-independent. 



C. Molecules for RENP 

Molecules are interesting candidates for RENP measurements from a variety of reasons: (1) 
they have a rich structure in energy levels, in particular energy level spacing of vibrational 
bands may give features of the photon energy spectrum helpful for identification of the 
RENP process; (2) closeness of different parity levels may enhance parity violating effects. 

In this Appendix we explain basic facts on molecules and present crude estimates related 
to molecular RENP, by taking I2 molecule as an example. In the fundamental Born- 
Oppenheimer approximation [17] the molecular wave function of diatomic molecules consists 
of three parts; electronic, vibrational, and rotational parts. Energy scale associated with each 
of these is clearly separated refiecting three different time scales of their motion. 

The vibrational part of wave function is usually derived after electronic wave functions of 
energy eigenstates are calculated taking the nuclear distance of two atoms fixed. For each 
fixed nuclear distance one has different energy functions corresponding to different electronic 
states. These energy functions are called potential curves. Potential curves have different 
equilibrium nuclear distances denoted by re (corresponding to different positions of energy 
minima) and different curvatures LOe at these minima, as illustrated in the conceptual diagram 
in Sec. 3. We need to compute matrix elements of electronic operators such as electronic spin 
S and electric dipole d. Even for two different electronic states these electronic operators have 
non-trivial overlap of vibrational wave functions because of different equilibrium distances 
and potential curvatures. This gives rise to the Pranck-Condon (PC) factor, as computed 
below. 

For simplicity we construct as a model of potential curves the Morse potential from molec- 
ular experimental data of Table 2 given in Sec. 3 and calculate the PC factors we need for 
RENP spectrum computations. 

C.l. Morse potential and vibrational energy eigenstates 

The Morse potential is a three-parameter fit to the potential curve and is given by 



with m the effective mass of two nuclei, in this case the half of I atom. This potential curve 
is plotted for X, A' and A states, three lowest electronic states in Fig.Cl, along with three 
vibrational levels. 




(CI) 




(C2) 
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Fig. CI The Morse potential for states, Fig. C2 Vibrational wave functions for 
X(v=15 indicated), A'(v=l) and A(v=2). Xv=15, A'v=l and Au = ~ 2 



FC product; AV-1 and Xv-15 



2 4 6 8 10 

Fig. C3 FC product between A' and X: the intermediate sum is taken up to Av (plotted 
abscissa is Av+1). 

This potential has analytic solutions for the energy eigenvalue problem: 

E.=Ue{v+^-)-^^{v+\f-De, V = 0,l,---[\-l], ^=^- (^3) 

The wave functions are given in terms of the associated Laguerre polynomial L". Examples 
of these wave functions are shown in Fig.C.l. 

C.2. FC factors 

The Franck-Condon (FC) factor is the overlap of vibrational wave functions, for instance 
between A'-A, and A-X. We need a product of FC factors between A'-A and A-X for RENP 
amplitude calculation. Result of FC product between A' and X for summation up to Av is 
shown in Fig.C3, which indicates that the approach to the final product value is fast, already 
seen Av ~ 2. 

C. 3. RENP spectrum 

The RENP spectral shape given in the text is used. We neglect the spin factor S expected 
to be order unity, and further ignore the presence of rotational levels. Examples of the 
photon spectrum function I{uj) are shown in Fig.C4 ~ Fig.C5. The molecular spin factor 
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12 A'v=l,2 NH and IH,ml=20mi:V 



12 A'v=l NH and IH: D vs M,ml=2()meV 




Fig. C4 I2 RENP spectrum between A' 
v=l and Xv=15. IH vs NH is compared, 
using the smallest neutrino mass 20 meV. 
The vertical scale is in arbitrary units. 
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Fig. C5 I2 RENP spectrum between 
A'v=l and Xv=15 in solid, and A' v=2 and 
Xv=14 in dashed. IH vs NH is compared, 
for two cases of the smallest neutrino mass 
2 (in dashed) and 20 (in solid) meV. The 
vertical scale is in arbitrary units. 



has been taken as a constant for simplicity. A more accurate calculation including the correct 
treatment of the spin factor is in progress [24] . 



C.4. The h A' state 

In this Subsection, we briefly describe how to produce the metastable A' states. I2 molecules 
form stable solids at room temperature with relatively high vapor pressure (~ 0.5 Torr). Fig. 
C6 shows saturated vapor pressure of iodine. [65] It is possible to prepare 1 atm I2 gas by 
heating up to about 450 K. Therefore, I2 is an easy-to-handle sample for both gas phase or 
matrix isolation experiments. Actually, many spectroscopic studies of I2 has been performed 
in both phases. In gas phase, however, at least three photons are needed to access the A state 
from X state. For example, X — )• ^ — )• a — )• A' may be possible by pure optical transition. 
As an alternative pathway, Koffend et al. used coUisional induced transition to prepare the 
A' state. [66] On the other hand, solid environment makes access to the A' state easy. It is 
known that A' — X transition is weakly dipole-allowed by interaction with a matrix. [67] 
This may be advantageous for preparing initial population in the A' state while both Ti and 
T2 in the solid phase are much shorter than in the gas phase. 



D. Coherence time measurements of para- hydrogen vibrational levels 

D.l. TRCARS method 

Solid para- hydrogen (PH2) is an attractive target to study the PSR phenomenon. Its prop- 
erties of high density (2.6 x 10^^ molecules cm~^) and long coherence time are well suited 
to observe explosive type PSR events, as discussed in Sec. 2. In this Appendix, we focus on 
the V = 2 (overtone) vibrational level of pH2, one possible initial state for PSR. Previous 
studies have reported that the v = 1 (of consisting H2 molecules) coherence decay exceeds 10 
ns [68-71], that is quite long as a condensed sample. A direct coherence time measurement 
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Fig. C6 Saturated vapor pressure of iodine. [65] 



is possible by investigating the temporal change in time-domain spectra. With this method, 
the coherence decay of the v = 1 state has been studied already [71]. Below, we describe 
our time-domain examination of the coherence decay of the overtone v = 2 state of solid 
para- hydrogen. 

TRCARS (time-resolvcd coherent anti- Stokes Raman scattering) [72] is a technique to 
trace the coherence decay in a real-time basis for Raman-type transitions. This is applica- 
ble to para-hydrogen vibrational levels since they are Raman-active. In this technique, the 
coherence produced by a short excitation pulse pair (pump and Stokes, whose energy differ- 
ence is equal to the Raman transition frequency ljq) are monitored by a delayed short probe 
pulse at t = At (sec Fig. Dl (a)). At t = 0, molecules in the excited sample are in a coher- 
ent superposition state between the upper and lower states of the Raman transition. This 
superposition state starts to decay after the excitation. At t = At, the probe pulse and the 
coherence that remains the sample generate the anti-Stokes pulse along the direction deter- 
mined by phase-matching condition. The intensity of this output pulse is proportional to the 
square of the coherence that remains at the delay At. In the case of usual exponential-like 
decays. 



where qo is the amplitude of the vibrational coordinate at t = 0, and r the coherence decay 
constant. The energy conservation law (Fig. Dl (a)) and the phase-matching condition (Fig. 
Dl (b)) are fulfilled in the TRCARS process as follows: 

^^L — ^^S = <^AS — Wp = Wo, 

kL - ks = kAs - kp = ko. 

Here and kj represent the frequency and the wavevector, related to each other by |ki| = 
n(LOi)cu}i, where n{ui) is the refractive index of solid para-hydrogen [73] at Wj. c the speed 
of light in vacuum. The subscripts L, S, AS, and P represent pump, Stokes, anti-Stokes, 
and probe, respectively. For v = 2, ooo is 8070.4(1) cm^^ in solid hydrogen [74]. ko is the 
wavevector generated at t = and then probed at t = At. 




(Dl) 
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Fig. Dl (a) TRCARS scheme for v = 2 coherence decay measurement in soHd para- 
hydrogen. The V = 2 coherence produced at t = by a pump (denoted by the subscript L) 
and Stokes (S) pulse pair is monitored by a probe (P) pulse delayed by t = At generating 
an anti-Stokes (AS) pulse, (b) The wavevector relations at t = (left) and t = At (right), 
(c) Experimental setup. Three nano-second pulsed lasers (Nd:YAG (wavelength: 355 nm) as 
pump, dye (498 nm) as Stokes, and another Nd:YAG (532 nm) as probe were irradiated to a 
sample of solid para-hydrogen. The phase-matched anti-Stokes signal (372 nm) was detected 
by a photomultiplier tube after spectral separation via a UV filter and a monochromator. 

The experimental setup used in this study is found elsewhere [75]. Briefly, the solid samples 
were obtained by the method described in Section 4.3.1. The sample temperature and the 
ortho hydrogen impurity concentration were controlled for optimization. For the TRCARS 
measurement as shown in Fig. Dl (b), three commercially available, nano-second pulsed 
lasers were used as the pump. Stokes, and probe sources. After spatial separation from the 
laser pulses, the anti-stokes signal pulse in the UV region was detected as a function of the 
delay. 

D.2. V = 2 coherence time of solid para-hydrogen 

Figure D2 shows the v = 2 coherence decay profiles under various conditions. Before opti- 
mization, the decay constants r were found to be 10-15 ns (Fig. D2 (b)-(d)). Ortho hydrogen 
impurities disturb the translational symmetry of the sample crystal, then an induced inho- 
mogeneous component contributes to the coherence decay [70] (Fig. D2 (b)). An additional 
contribution occurs from elastic scattering of the vibrationally excited molecules by ther- 
mal phonons [69, 76] (Fig. D2 (c)). This contribution is known to obey a law on the 
temperature dependence of the decay rates [77] . We observed this dependence for the v = 2 
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Fig. D2 Coherence decay profiles of the vibrational states v = 2 of solid para-hydrogen, 
(a) The decay measured under the optimized condition (after annealing, at the sample 
temperature T = 4 K, and at the ortho hydrogen concentration o/p = 0.01%). The others 
were obtained under the same condition except for (b) o/p = 1.9%, (c) T = 8 K, and (d) 
before annealing. Each profile is normalized by its peak value. Solid lines are exponential 
fits to the data, in which the data points corresponding to the normalized intensity below 
0.1 were used. The decay constant r of (a) is found to be ~ 50 ns. 



excited state. [75] Finally, the intrinsic structural defects induced during sample preparation 
accelerates the coherence decay (Fig. D2 (d)). 

By improving the sample conditions described above, we obtained quite a long decay 
profile (Fig. D2 (a)). Here, the ortho concentration was lowered to 0.01 %, and the sample 
temperature was cooled down to T = 4 K. To optimize the decay further, an annealing 
process for half an hour at T = 9 K was performed before the measurement. As seen in Fig. 
D2 (a), the coherence was kept almost constant up to ~ 50 ns just after the excitation, . This 
can be accounted for by a transient stimulated Raman process induced on the excitation. 
[78] The profile after this plateau has an exponential- like decay. Fitting this part to Eq. Dl 
gives r ~ 50 ns. A deviation from an exponential function at longer delay (also seen in Fig. 
D2 (d)) is discussed elsewhere. [75] The obtained coherence time for v = 2 is comparable to 
that for = 1. In the PSR process, a coherent superposition state of the upper and lower 
states develops in time. It is required that the disturbance of the coherence is weak during 
the process. The long coherence times of v = 1 and v = 2 show that a nano-second scale 
measurement can trace PSR signals for solid para-hydrogen when these states are employed. 

Generally in solid samples at low temperatures, the dephasing mechanism in the gas 
phase as Doppler broadening and collisional relaxation is minimized. Instead, the inter- 
atomic/intermolecular distance is close as a few A, and the inhomogeneity and the 
fluctuation of the solid crystal structure causes corresponding variations in the inter- 
atomic/intermolecular interaction. These variations result in fast relaxation usually in a 
time scale of 10-100 ps. However, the vibrational states in solid para-hydrogen have such a 
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long coherence time as a few 10 ns, because of a rather weak intcrmolecular interaction and 
the high homogeneity due to large zero-point vibration around lattice points. 

The coherence time of 50 ns obtained here for the v = 2 state of solid para-hydrogen sat- 
isfies the assumption in a simulative work of PSR. [14] The coherence time can be further 
improved by a single crystal sample by the pressurized liquid method [79], that was exhib- 
ited for V = l.[38]. To succeed in PSR observation, the excitation sample density is another 
important factor. Additionally, the initial coherence preparation by intense lasers are nec- 
essary to observe an explosive PSR. In the present solid para- hydrogen experiment, the low 
damage threshold of our sample crystals (~ 1 MW/cm^) limits the maximum laser inten- 
sity. Crystals with high damage threshold can be obtained by the pressurized liquid method 
mentioned above. In order to achieve the initial coherence of 0.1 — 0.01 as achieved in the 
V = 1 excitation, [80] the excitation lasers are required to have a good coherence. These laser 
system, as discussed in Sec. 4.2.2, is being exploited for PSR in gas para-hydrogen experi- 
ment, which is a promising candidate for explosive PSR because the samples are basically 
damage-free. 



E. Experimental studies on PSR/RENP targets in condensed phases 

Atomic and molecular systems, examined in our group for PSR/RENP targets other than 
pll2, I2, and Xe, are briefly introduced. Condensed matter targets are advantageous in the 

following three aspects. 1) High density energy storage is possible in a typical number den- 
sity of 10^^ cm^^ for target atoms and molecules. 2) Coherence control is possible by tuning 
inter-atomic or inter-molecular interaction to suppress undesired processes for decoherence 
accompanied by dissipation of energy. 3) Doppler broadening in the spectral linewidth is min- 
imized by keeping the solid sample at low temperatures. Despite a disadvantage for some 
materials of relatively low damage threshold to laser exposure, solid targets are promising 
for future experimental research for PSR/RENP. In this Appendix, spectroscopic character- 
ization of condensed-phase atomic and molecular systems are described, i.e., HE molecules 
in solid pH2, atomic N in a carbon cage of Ceo, and Bi atoms in solid Ne. 

E.l. Bismuth in neon matrix 

E.1.1. Energy levels of Bi. As a heavy element with large spin-orbit coupling, bismuth 
is one of the excitation targets for the observation of radiative emission of neutrino pair 
(RENP). Figure El shows the energy levels for Bi atom. Similarly to the other group-15 
elements, the low- lying excited states, ^P3/2,i/2 a-^id ^P5/2,3/2; a-i'e stemming from the same 
electron configuration, 6s^6p^, as that for the ground state, ^83/2- Due to the large spin- 
orbit coupling, however, unlike the lightest group-15 element of nitrogen, different J states 
for bismuth are largely separated, i.e., 10504 cm~^ for ^P (J=3/2 and 1/2) and 4019 cm~^ 
for (J=5/2 and 3/2). As a result of the mixing of different spin states, orbital angular 
momentum, L, and spin angular momentum, 5, are no longer good quantum numbers. As a 
result of the mixing of different spin states by the spin-orbit interaction, transitions become 
possible for a pair of states with the same electron configuration. The radiative lifetimes for 
these transitions, ~10~^— 10~^ seconds in the gas phase, are long enough to be detected by 
using an instantaneous excitation with a nanosecond laser pulse. These low-lying excited 
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Fig. El Level diagram of bismutli atom (gas pliase). 



states were observed in solid Ne matrix [81]. Using higher excited levels accessible by non- 
resonant two-photon excitation, single-photon super-radiance was observed with Bi atoms 
in the gas phase [82]. 

E.1.2. PSR/RENP transition. In Fig. El, two energy levels are noted. One is a metastable 
^05/2 state (6p^) and the other is the ground ^83/2 state (6p^). These two states can be 
connected via intermediate states at higher energies, e.g., ^P5/2 (6p^7s^), with a combination 
of El and Ml (practically El) transitions. The set of these energy levels is an example for 
a A-type ladder for paired super-radiance (PSR) using bismuth atoms. If the coherence is 
developed between the ground state and the metastable excited state, and if an appropriate 
trigger source is applied to the system, explosive two-photon emission of radiation by PSR 
or even RENP can be promoted. A PSR rate of ~17 kHz was calculated for the number 
density 10^® cm~^ of excited bismuth atoms. 

E.1.3. Matrix isolation technique. In order to realize PSR or RENP, to produce a system 
containing a large number of excited atoms in a small volume is crucial. The use of a solid- 
state material is demanded for realization of a large atomic density, thus for prolonged 
observation in the detection of events. Matrix isolation technique has been a spectroscopic 
tool applicable to reactive species such as atoms, molecules, open-shell radicals, and clusters 
of metal atoms. The target species are entrapped in small vacancies in rare-gas solids of Ne, 
Ar, Kr, or Xe. The number density of target atoms can be increased up to lO^'' cm~^ under 
well-isolated conditions. For matrix-isolated species, relaxation pathways to the long-lived 
metastable state can be enhanced, by which the excited-state population is increased to a 
substantial fraction or even inverted against the ground state, as recently demonstrated for 
Yb atoms in solid Ne [83]. 

E.1.4. Laser induced optical emission. We performed leaser-induced optical emission 
spectroscopy to determine energy levels and lifetimes for Bi atoms in solid Ne matrix. Vapor 
of bismuth emanated from a bismuth-containing molybdenum crucible at ~1000 K was co- 
deposited on a cold surface of sapphire at 3 K together with an excess of neon gas. After 
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Fig. E2 Emission spectra of Bi atoms trapped in a solid Ne matrix at 3 K using the 
excitation with nanosecond laser pulses at 300.7 nm (upper trace) and 302.0 nm (lower 
trace) for the ^P3/2'^^S3/2 transition. 

the deposition, nanosecond laser pulses tuned at transition wavelengths for Bi atom were 
irradiated on the solid sample in a grazing angle. The emitted light was dispersed by using a 
grating spectrometer (Acton SP300i) and detected by using a CCD camera (PI SPECIO) for 
the observation of emission spectra. A photomultiplier (Hamamatsu R928) was used for the 
measurement of lifetimes. In order to reduce stray light of the excitation laser, a long-pass 
optical filter (Schott Glass Filter) was used. 

First we excited the ^P3/2'^^S3/2 transition. Figure E2 shows typical emission spectra for 
the excitation at 300.7 nm. Emission bands corresponding to three relaxation pathways were 
observed. Two of them were the transitions from the upper state of ^P3/2 to the metastable 
states of ^D5/2 (565 nm) and ^03/2 (460 nm). The other one was the transition from the 
metastable state of ^05/2 (645 nm) to the ground state of ^83/2- Each band was composed 
of a few lines, for which the relative intensity varied upon different excitation wavelengths. 
At the excitation wavelength of 300.7 nm providing maximum emission intensity, relatively 
sharp emission lines were intensified at the blue edge of the emission band. A broad emission 
band at 570 nm in the lower trace in Fig. E2 may be associated with another transition 
of ''Pi/2^^D5/2 (6p27s^^6p3), whose upper sate of ^Pi/2 locates slightly below the ^P3/2 
(6p^) level. The lifetime for this transition at 570 nm was several times shorter than that 
for the transition at 565 nm. 

Figure E3 depicts decay profiles for the three transitions for Bi atom in solid Ne. The single 
exponential fits provide lifetimes of 1.05 ms for ^05/2— >-^S3/2, 0.31 ms for ^P3/2— ^^05/21 
0.35 ms for ^P3/2~^^D3/2- These lifetimes in a Ne matrix are two orders of magnitude shorter 
than those in the gas phase. The latter two, corresponding to relaxation to the metastable 
^^5/2, 3/2 states, showed a single exponential decay with a characteristic lifetime of ~0.3 ms, 
while the former one, form the metastable state to the ground state, showed a characteristic 
feature of rise in ~0.2 ms and a decay in a longer lifetime of ~1 ms. These profiles are 
consistent with a picture that the population in the metastable state, ^05/2, increases via 
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Fig. E3 Decay of emission bands for Bi in solid Ne witli the excitation at 300.6 nm, 
205/2^^83/2 (645 nm), ^P3/2^^B,/2 (565 nm), and ^P3/2^^By2 (460 nm). 
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Fig. E4 Excitation profiles for Bi in solid Ne, a) 2p3/2^^S3/2 (~300 nm) and b) 
^^1/2^^83/2 (~460 nm). Inset shows a typical emission spectrum with fitting curves by 
three components. 

the radiative decay of ^P3/2~^^D5/2 in a time scale of ~0.3 ms and decreases via that of 
^D5/2-^'^S3/2 in ~1 ms. 

In order to locate the energy levels of 2P3/2 and states, excitation profiles were 

examined. By using the photomultiplier detector equipped to the spectrometer, emission 
intensity at a fixed wavelength was measured and plotted as a function of laser-excitation 
wavelength. Figure E4 shows the excitation profiles for Bi atoms in solid Ne. Corresponding 
to the absorption for 2p3/2-(— ^83/2, the maximum was found at 300.5±0.25 nm by plotting 
the emission intensity at 565 nm as shown by closed circles in green in Fig. E4(a). The 
bandwidth for this transition was ~2 nm (~220 cm~^). The excitation profile for the 570- 
nm emission band (see the lower trace in Fig. E2) was much broader as shown by open circles 
in Fig. E4(a). For 2Pi/2-^^S3/2, the maximum was found at 458.5ib0.25 nm by plotting the 
emission intensity at 645 nm as shown by the three profiles in Fig. E4b, for which spectral 
decomposition was performed for each of the observed emission spectra at different excitation 
wavelengths and the excitation profile was plotted for each of the three components in the 
emission spectrum. The bandwidth for the ^Pi/2'^^S3/2 transition was ~4 nm (~190 cm"^). 
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Table El Observed transitions for Bi atoms in solid Nc matrix and comparison with the 
transition energy in the gas phase. "Radiative hfetimes in sohd Ne with an excitation at 
300.6 nm. ''Difference of term energies. 



Transition 



Wavelength/nm Lifetime/ms" Energy/cm 



-1 



Energy(gas) / cm' 



-1 




300.5 - 33280 

458.5 - 21810 

459.5 1.05 21760 

564.7 0.31 17710 

643.2 0.35 15550 



33165 

21661 

21746'' 

17727^* 

15438 



E.1.5. Linewidth and coherence. At the blue edge for each of the emission bands in Fig. 
E2, a relatively sharp peak is discernible. This peak corresponds to the zero-phonon line for 
each transition. The linewidth for these relatively sharp lines was ^^10 cm~^, while the total 
width for the bands was ~200 cm^^. Even for the narrow lines, the linewidth corresponds 
to a natural lifetime on the order of ^lO^^'^ seconds. Since the observed lifetimes in Fig. 
E3 is much longer, ~1 ms, the observed spectral linewidth is ascribed to inhomogeneous 
broadening due to interaction of the target atoms with the matrix media. With substantial 
deviations in the transition wavelength, it is considered to be difficult to develop coherence 
longer than a picosecond in this system. 

The laser-induced emission spectroscopy revealed energy levels for all the 6p'^ states of Bi 
atom in solid Ne. The observed transitions in this work are summarized in Table El. Term 
energies for the four low- lying excited states of Bi atom in solid Ne remain intact as those in 
the gas phase within a bandwidth in the matrix spectra. Lifetimes were revealed to be two 
orders of magnitude shorter than those in the gas phase. Despite the advantage of the long 
lifetimes of ~1 ms, the apparent inhomogeneous broadening would suppress coherence time 
shorter than expected. To remove undesired interaction between embedded guest atoms and 
hosting matrix atoms is crucial for realizing a system for PSR or RENP. 

E.2. HF molecule trapped in solid pH2 

E.2.1. Matrix isolation spectroscopy using solid pH2. Spectroscopy of molecules embed- 
ded in condensed phases is called matrix isolation spectroscopy, which has long been used to 
pilot gas phase spectroscopy. Not only unstable but also stable molecules have been subjected 
to studies by this method. Rare gas matrices such as solid Ne and solid Ar have been widely 
used because of their chemical inertness and of relatively weak perturbative interactions. 
The perturbation in solid rare-gas matrices sometimes makes spectral linewidths broader 
compared to those in the gas phase spectroscopy, to conceal detailed spectroscopic informa- 
tion. Small molecules in solid pH2, on the other hand, exhibit extremely sharp lines to reflect 
quantized rovibrational states of the entrapped molecules in the condensed phase. The full 
width at half maximum (FWHM) for the vibrational 1/4 band of methane, CH4, in solid pH2 
is as narrow as 0.015 cm~^, whose signal is more than one order of magnitude sharper than 
that observed in conventional solid rare-gas matrices [84]. Such a narrow linewidth of ~0.01 
cm~^ is a promise of a long relaxation time advantageous for the observation of coherent 
phenomena. 
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Fig. E5 IR absorption spectra of HF molecules in solid pH2 for the sample a) as prepared 
and b) 1-day after preparation. 



E.2.2. Preparation of chemically doped solid pH2. To prepare the doped solid pH2, the 
closed cell method (Sec. 4.2.1) can be used. Premixed pH2 gas containing small amount of 
the target molecules is introduced into the cooled closed cell and a transparent doped crystal 
grows in the same manner as pure pH2. However, only dopants of small intermolecular inter- 
action like methane are isolated by this method while molecules with strong interaction like 
HF are difficult to be doped dispersedly. In contrast, the rapid deposition method developed 
by Fajardo et al [85] can isolate a number of molecular species in solid pH2. They include 
CH4, CH3F, CH3I, CH2CII, CO, CO2, NO, N2O, H2O, HF, and HCN. Several millimeters 
thick, optically transparent solid pH2 can be prepared by this method. The pH2 and dopant 
molecules are co-deposited onto a cold substrate (typically BaF2 at 2~4 K) placed in a 
cryostat. The substrate is contacted on the copper block with an indium gasket. The crystal 
grows perpendicularly from the substrate and a polycrystalline aggregate of h.c.p. and f.c.c. 
having their c axes roughly normal to the substrate is produced. To remove the metastable 
f.c.c. crystal, the mixed crystals are heated up to ~ 5 K for 10 to 30 min. This annealing 
aligns c axes normal to the substrate. 

E.2.3. Infrared spectroscopy of HF in solid pH2- Hydrogen fluoride, HF, in solid pH2 has 
been studied in our laboratory as a target for super-radiance in the condensed phase [86]. 
Super-radiance was observed indeed in the gas phase for the rotational transition of HF, 
where gaseous HF was optically pumped into the rovibrationally excited state {v=l, J) 
for realizing total population inversion against the lower state {v =1, J-1), resulting in the 
super-radiance with far infrared (FIR) emission of photons [21]. 

We measured IR absorption spectra for HF molecules in solid pH2 at 3.6 K by using an 
FTIR spectrometer. The observed spectra for HF/pH2 are shown in Fig. E5. Based on the 
systematic spectral change caused by the difference in HF/pH2 concentration (not shown), 
we assigned the observed IR absorption lines in the region of 4000-3000 cm~^ to isolated 
molecular HF, dimers (HF)2, trimers (HF)3, and clusters (HF)„ (n>4) as indicated in Fig. 
E5. 

Time evolution in the IR absorption spectra in Fig. E5(a) and E5(b) indicated that the HF 
molecules could migrate in the solid pH2 at 3.6 K, as was noted by the decrease in the line 
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Fig. E6 Level diagram of nitrogen atom (gas phase). 



at 3970 cm~ for isolated HF and the concomitant increase in the bands at lower frequencies 
for clusters (HF)„. The aggregation of HF molecules in solid pH2 indicates that the isolated 
HF molecules cannot be the targets for long-term observation. Moreover, the linewidth for 
the isolated HF in solid pH2 is rather broad, ~4 cm~^, probably due to matching between 
rotational energy for HF and phonon energy for solid pH2. Note that the linewidth for the 
signal at 3820 cm^^ for the dimer, (HF)2, is much narrower than that for the isolated HF, 
probably because the rotational motion is quenched. Solid pH2 is, in general, a good matrix 
medium for trapping molecules with minimal interaction, whereas, for the small but polar 
molecule such as HF, another difficulty may arise in migration or rotation-phonon coupling. 



E. 3. Nitrogen atom in fullerene Cqq 

E.3.1. N@Cqo: discovery and characteristics. Fullerene Ceo is a hollow, closed-cage 
molecule of carbon with a dimension of ~1 nm, in which atoms of the other elements can be 
accommodated. Atoms of lanthanides and some transition metals are encapsulated in larger 
fullerenes such as C82 by carbon arc using a metal-carbon composite rod as an electrode. 
Encapsulation of atoms in a Cgo cage of icosahedral symmetry has turned out to require 
another method, i.e., ion implantation. Li@C6o [87] is a prototype of ion- implanted fullerene 
molecules (MQCeo depicts M inside a Ceo fullerene cage) [88], while NQCeo is an unique 
member among the endohedral fullerenes [89]. 

Unlike metallofuUerenes having a metal atom locating close to inner walls of the carbon 
cage, the nitrogen atom in Ceo is believed to be located at the center of symmetry of the 
Ceo molecule. Moreover, the electron spin, S = 3/2, of atomic nitrogen is retained even in 
the carbon cage. Furthermore, the phase coherence time in the electronic ground state can 
be longer than 0.2 ms as revealed by electron paramagnetic resonance (EPR) [90]. With all 
these characteristics in mind, one can believe that the interaction between the trapped atom 
and the hosting cage should be minimized, thus enjoyed by further experimental studies. 
Once a crystalline form of pure NQCeo is obtained, the density of atomic nitrogen amounts 
to 2.0x10^^ cm~^, each atom being well isolated as a spin carrier for coherent phenomena. 
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E.3.2. Electronic states of atomic nitrogen. Figure E6 illustrates atomic energy levels for 
some low-lying electronic states of free nitrogen atom. The excited states, ^P3/2,i/2 and 
^1^5/2,3/2) as well as the ground state, ^83/2, are stemming from the same electron con- 
figuration of 2s^2p'^. Therefore, transition between any pair of states is forbidden by El 
mechanism. To excite N atom in its or states, some higher excited states to which 
the El transition is allowed from the ground state should be involved. Note that transitions 
between the low-lying electronic states are possible by Ml and/or E2 mechanism, though 
the rate of these transitions is relatively slow. 

The excitation channel by El mechanism is realized by vacuum-UV transitions of 3s-<— 2p 
(120 nm and 113 nm). Pumping to the higher level of ^P5/2,3/2,i/2 is followed by relaxation to 
the low-lying excited state of ^P3/2,i/2 or ^05/2,3/2 via transitions induced by spin-orbit inter- 
action. When the incoherent energy loss is minimized, population at the low-lying excited 
state increases by repetitive pumping at a high rate. 

E. 3. 3. Preparation of N@Cqq . Formation of NQCeo is conducted under vacuum where Ceo 
molecules are bombarded with positively charged ions of nitrogen, which are accelerated to 
have a kinetic energy of several tens of electron volts. Eventually, some are trapped inside the 
Ceo cage and extracted by solvents as a neutral molecule of N@Ceo- The raw material after 
ion implantation contains at most one thousandth of N@Ceo in Ceo- Therefore, separation 
of N@Ceo from Ceo is crucial. 

We succeeded in isolation of NQCeo using a recycling-HPLC system. Figure E7 shows a 
chromatogram for the final step of separation of NOCeo and Ceo- Due to the faint difference in 
interaction between the molecule and the surface of the porous material in the 5PBB column, 
it takes a few percent longer times for N@C6o to pass through the column than those for 
Ceo- After the repetitive passage of four cycles, the fraction of N@Ceo was collected for the 
measurement of optical properties. 

The inset in Fig. E7 shows UV-vis absorption spectrum of N@Ceo after purification. The 
three major bands in the UV are essentially the same as those for empty Ceo, indicating 
that the interaction between the N atom and the Cqo molecule is not discernible at all. Since 
the excited states for NQCeo, ^P and ^D, are not accessible directly by the transition of El 
mechanism, corresponding absorption lines are not observed. It is natural to consider that 
the low-lying excited states of atomic nitrogen in a Ceo cage remain intact as those in the 
gas phase. 

E.3.4- EPR detection. Here, the EPR properties in the ground state are described briefly 
to demonstrate the significance of N@Ceo- In the inset in Fig. E8, three narrow lines by 
cw-EPR measurement correspond to the transitions between different hyperfine states in 
^83/2 of ^^NQCeo at room temperature in a solution of carbon disulfide. Among the total 9 
lines for the EPR transitions, three lines of constant M/ = +1,0, or —1 arc overlapping in 
each of the three lines in the spectrum. This indicates that there is no noticeable anisotropy 
for the electron spin in NQCeo- 

The main panel in Fig. E8 displays free induction decay (FID) for the pulse-EPR signal 
for a powder sample of N@Ceo/C6o mixture, tuned close to the transition of the M/ = -|-1 
line at lower field. The moderate wiggle of a period of ~3 jis is due to detuning from the line 
of Ml = -1-1, while the rapid modulation is due to interference by the signals of Mj = and 
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Time /us 



Fig. E7 HPLC chromatogram for sepa- 
ration of N@C6o and Ceo, i-e-, the 6th recy- 
cling step using 5PBB {(f)10x250 mmx3 
columns), toluene (3.0 mL min~^), and a 
UV detector (333 nm). Inset shows UV- 
vis absorption spectrum of the fraction at 
2h52m in the main panel. Using absorption 
coefficient for Ceo (4.8x10^ L mol^^ cm~^ 
at 330 nm) and the volume of the solution 
(2.5 mL), the amount of isolated NQCeo is 
estimated to be 3x10^^ molecules. 



Fig. E8 FID signal for the Mj = +1 
transition of ^^NOCeo (/ = 1) obtained 
from X-band pulse-EPR measurement at 
room temperature for the powder sample 
of N@Ceo/C6o mixture (main panel). Inset 
shows X-band cw-EPR spectrum of purified 
N@Ceo in carbon disulfide. The triplet with 
a hyperfinc splitting of 15.8 MHz belongs 
to ^"^NQCeo (/ = 1), while the weak dou- 
blet with a splitting of 22.1 MHz belongs to 
i5N@Ceo {I = 1/2). 



Mi = —1 which are separated by 15.8 and 31.6 MHz from the Mi = +1 line, respectively. 
Owing to the narrow linewidth even in a solid state, coherence can be maintained at least 
for several microseconds. 

In conclusion, we have developed a system for preparation of purified N@Ceo in a micro- 
gram order, ~10^^ molecules. UV-vis absorption spectra indicate that the nitrogen atom in 
a Ceo cage stays intact as that in the gas phase. The excited-state properties are intriguing 
in view of the coherent phenomena induced by radiation field. 
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